Chapter 2 

The Role of Radioactive Isotopes in Astrophysics 

Donald D. Clayton^ 



2.1 History of Nucleosynthesis and Radioactivity 

Radioactive nuclei can not be not infinitely old or they would now be gone. This 



point had first been made by Lord Rutherford, who concluded in 1929 (Rutherford 



1929) 1 that the elements had been created in the sun 100 Myr ago and had somehow 



got from sun to earth! Terrible astrophysics, but his argument showed profound 
appreciation of the conundrum of live radioactivity on earth; namely, the nuclei must 
be created. By the middle of the 20 th century, the consequences had progressed no 
further. The question, "Why aren't they all gone?" went largely unspoken. 

Radioactivity had been discovered by Henri Becquerel in 1896 when he showed 
that the element uranium emitted radiation that would develop images of the ura- 
nium salts on an adjacent photographic plate. Marie Curie coined the term radioac- 
tivity to describe this phenomenon after she demonstrated that other elements pos- 
sessed the same spontaneous property. Rutherford separated the rays into three types 
- a, j3, and 7- by their physical properties. These are described in Chapter 1. 

In mid-20 th century, astronomers discovered metal-poor stars in increasing num- 
bers, especially after the introduction of the CCD chip, which greatly increased the 
sensitivity of spectroscopic observations of stars. The abundances relative to hy- 
drogen in old low-mass stars reveals the initial composition of the gas from which 
those stars formed. Observed metallicities ranged over a continuum of values from 
(Fe/H)/(Fe/H)© = 10~ 4 to 1 and even greater, so the synthesis of metals must have 
occurred after the Big Bang. Today we know that creation of the chemical elements 
occurs in stars by synthesis of them from initial H and He. By the same token, 
natural radioactivity must be the product of nucleosynthesis in stars. 
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2 2 The Role of Radioactive Isotopes in Astrophysics 

2.1.1 Two very different Pioneers 



There could hardly be two more different scientists than the two that made the first 
significant steps toward nucleosynthesis in stars of the elements heavier than carbon. 
Both were Europeans disrupted in careers by world war II. The originator of nucle- 
osynthesis in stars arose from non-privileged Yorkshire mill villages. The pioneer 
in glimpsing how nuclear magic numbers can be observed in isotopic abundances 
arose from a distinguished Viennese family of geologists. 



Sir Fred Hoyle (1946, 1954) 



In a flash of astrophysical brilliance, Fred Hoyle created the theory of nucleosyn- 
thesis in stars with two monumental papers (Hoyle et al. 1946 1954) . The first 



( |Hoyle et al.| Q946) demonstrated that stars will naturally evolve to central tem- 
peratures 4x 10 y K and greater, and that nuclear statistical equilibrium (NSE) at 
high temperature T creates an abundance peak at 56 Fe that Hoyle interpreted as 
the explanation of that observed abundance peak. This explains why a graph of the 
abundances of chemical elements in the sun as a function of atomic weight declines 
almost monotonically from H to Sc (A=l to 45), whereafter it increases sharply to a 
broad abundance peak centered on A=56 (Fe). Readers unfamiliar with the meaning 
and issues of abundances in the universe can find an eloquent essay in the introduc- 
tory chapter of Clayton ( 2003)>. 

Hoyle's second paper (Hoyle et al. 1954} described how the elements from car- 



bon to nickel are synthesized during the advanced evolution of massive stars. This 
is called primary nucleosynthesis because the elements are fused from the initial H 
and He of the massive stars. Hoyle described how the ashes of each thermonuclear 
burning epoch became, upon contraction and heating, the fuel for the next group of 
elements synthesized. When its core collapses, much of the overlying new elements 
is returned to the interstellar medium. Such nucleosynthesis creates new heavy nu- 
clei and is responsible for increasing the metallicity of the interstellar medium as it 
ages. 

Because the origin of our chemical elements is one of the grand theories of sci- 



ence, ( Clayton 2007 (went to some lengths to create what he called Hoyle 's equation 



as determined from careful reading of his 1954 paper (Hoyle et al. 1954} . Hoyle's 
words and calculations pointed clearly to ideas of nucleosynthesis in stars that he 
was advancing for the first time and that are more sweeping than detail-oriented se- 
quels. Hoyle's discussion was phrased in terms of the mass Am new of new primary 
isotopes that are being ejected from massive stars. His basic approach to stellar nu- 
cleosynthesis is 



dm„ 



dt 



— Hnucl 



where 



H nuc i =B{t')E\(t' -t)Z k Am k 



(2.1) 
(2.2) 
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is the Hoyle nucleosynthesis rate. B(f') is the stellar birthrate of stars having total 
mass such that they evolve to end their lives at time t, Ev(f',f) is an operator (rather 
than a number) that expresses the nuclear and stellar evolution during its lifetime 
from t' to t, and Em^ is the mass of isotope k ejected at time t. Then a sum over 
all presolar birthdates t' selects the appropriate stellar masses. Stellar evolution was 
only dimly perceived in 1953 when Hoyle wrote this paper. The structure of red 
giants was its current literature frontier, and Hoyle was perhaps the world's leading 
expert on its ideas, enabling him to discern the more fruitful subsequent evolution 
that he advanced in this 1954 paper. 

Hoyle subtitled this paper The synthesis of elements from Carbon to Nickel. These 
primary isotopes are ejected from massive stars that progress through a series of 
core evolutions that Hoyle laid out for the first time. Hoyle explained that gravita- 
tional contraction causes temperature increases after each nuclear fuel is consumed, 
and he described the nuclear burning during each advanced core evolution. Because 
those massive stars all evolve almost instantaneously in comparison with galactic 
timescale, Hoyle takes Bu> (t ) to be the birthrate of all such massive stars at time t , 
and it clearly equals their death rate at the same time if the numbers of stars are to 
change only slowly. The subscript M > characterizes stars too massive to become 
white dwarfs. For those large-mass stars Hoyle et al. (19541 predicted that collapse 



of the final central evolved core is inevitable. For those massive stars Hoyle's equa- 
tion expresses the rate of ejection of new primary isotopes from C to Ni as 

dm(C-Ni) 

y — '- = B M> (t)Ev nuc ,Z k Am k (2.3) 

at 

For these short-lived stars, Ev is an operator (rather than a number) that expresses 
the nuclear and stellar evolution during the stellar lifetime. It replaces the initial 
composition of the star by the composition it has attained at the time when its core 
collapses. Hoyle attributed the mass and identity k of new primary isotopes ejected 
per massive star to the following successive core burning phases: 12 C and lfi O from 
He burning; 20 Ne, 23 Na, and 24 Mg from C burning; additional lfi O and 24 Mg from 
Ne burning; 28 Si and 32 S from O burning; 32 S, 36 Ar and 40 Ca from photoalpha 
reactions on 32 S and heavier alpha nuclei during later heating of O-burned matter 
by the inevitable contraction; and finally 52 Cr, 56 Fe, 60 Ni from subsequent nuclear 
statistical equilibrium. Hoyle also correctly stated that neutrino emission governs 
the collapse timescale when core temperature exceeds T=3x 10 9 K. 

Hoyle's equation expresses a breathtakingly modern view of the metallicity- 
increasing nucleosynthesis during galactic history. Hoyle missed only the complete 
photonuclear rearrangement during Si burning and the n/p ratio within the NSE. 
But his equation, given above, remains correct today. Like Schroedinger's equation, 
for which much work is required to determine the appropriate Hamiltonian opera- 
tor used within it, so Hoyle's equation involves much work to determine the Hoyle 
nuclear evolution H nuc i. 

Hoyle et al. ( 1954] > also distinguished the idea of secondary nuclei, those whose 



created abundance derives from initial seed concentrations of primary nuclei that 
had been created in previous massive stars, thereby seeding the interstellar matter 



4 



2 The Role of Radioactive Isotopes in Astrophysics 



with primary elements. These seed nuclei are required in order that secondary nu- 
clei be produced from them in stars. He emphasized especially 14 N, ls O, 19 F and 
22 Ne in that regard, each of which depends on initial primary C and O nuclei within 
the initial composition of later-generation massive stars. Their yields of secondary 
nuclei do not obey Hoyle's equation but instead are proportional to the initial metal- 
licity of each star. Hoyle's paper also first noted that 22 Ne would be a source of free 
neutrons; indeed, it is today their major source in burning shells of massive stars, 
although that insight is usually attributed to later emphasis by others. 

It is unfortunate that Hoyle did not put to paper the equation he envisioned and 
described. Had he done so, clearer scientific visibility of his unparalleled achieve- 
ment would have followed more easily. Research in nucleosynthesis has followed 
his approach during five decades of countless numerical models evaluating Hoyle's 
equation. Nonetheless, the astrophysical world began to inappropriately cite a paper 
to appear, with Hoyle as coauthor, in 1957 for the general origin of nucleosynthesis 
in stars rather than Hoyle's works. Hoyle's great achievement lay somewhat ob- 
scurred until modern times, at least in terms of citations of research works. 

Hans Suess, Nuclear Shell Model, and Abundances 

Hans Suess was born in Vienna in 1909 to a distinguished family of geologists. His 
chemical training and focus on abundances of the elements in meteorites enabled 
him to establish Suess 's rules for the abundances and their isotopes. The elemental 
abundances were not well known after WWII, but testing of a theory of nucleosyn- 
thesis needed that data base. In the late 1940s Suess began to argue that certain 
regularities of the abundances had to reflect nuclear properties of their isotopes. He 
used those systematics to coauthor with Hans Jensen a discovery paper ( |Haxel et aL 
1949 ) for the magic numbers of nuclear shell structure. These are hugely important 
for nucleosynthesis. The magic numbers are favored numbers of identical nucleons, 
either protons or neutrons, in the sense that clusters of those numbers of identical nu- 
cleons have larger than normal binding energy. They are 2, 8, 14, 20, 28, 50, 82 and 
126, and they result from the combined effect of a deep spherical potential binding 
the nucleons together, coupled with a very strong spin-orbit energy that moderates 
energy gaps between differing nuclear shells. For that paper Jensen shared the 1963 
Nobel Prize in Physics with Maria Mayer for their independent theoretical work on 
how the spin-orbit force establishes the magic numbers. Suess, working in Hamburg 
and Heidelberg, had seen the evidence of magic numbers in the abundance regular- 
ities shown by nuclei having a magic structure of neutrons or protons. For example, 
the abundances of the isotopes of barium implicated neutron irradiation in their nu- 
cleosynthesis. Suess amazingly divined that the very high abundance of 56 Fe had to 
reflect the doubly-magic properties of its 56 Ni isobar having A^=28 for both protons 
and neutrons. That correct assertion was not accepted for a decade because Hoyle's 
papers (and also B 2 FH) had maintained that that abundance peak was established 
within a nuclear equilibrium having excess of neutrons relative to protons, so that 
the abundance of Fe isotopes could reflect their own nuclear properties rather than 
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those of Ni. This error persisted for a decade, and had many dead-end astrophysi- 
cal consequences. These have been intensively reviewed ( |Clayton 1999), who also 
reviews Suess's pioneering papers from the 1940s. 

The magic numbers were also to play a pivotal role in the theories for nucleosyn- 
thesis of the elements heavier than the Fe peak. Some of the earliest testaments to 
the correctness of the theoretical ideas lay in their interplay with the magic num- 
bers. Suess himself called for two neutron-addition processes needed to account for 
magic-number abundances within the heavy elements. After moving to the United 
States in 1950, Suess coauthored with Harold Urey an immensely influential tab- 
ulation ( |Suess and Urey 1956| l of the abundances of the elements, largely from 
chemical analyses of meteorites rather than from astronomical observations. That 
review paper by Suess & Urey became a cornerstone of empirical evidence for nu- 
cleosynthesis in stars. 



2.1.2 The Second Decade 

Hoyle's two papers were not enough to establish the theory of nucleosynthesis in 
stars. Contributions from other pioneers came together to launch the full theory and 
to win the acceptance that it enjoys today. 



William A. Fowler and Burbidge, Burbidge, Fowler and Hoyle (B 2 FH) 

William A. Fowler had met Fred Hoyle in 1953 when Hoyle famously predicted the 
existence and the energy of the + second excited state of 12 C. In that first applica- 
tion to nuclear physics of what has come to be called the anthropic principle, Hoyle 
had argued that if such a nuclear state of 12 C did not exist, neither would we! Fowler 
said that that prediction was what "really hooked me on nucleosynthesis". Fowler 
met Hoyle again during his 1955 sabbatical leave in Cambridge UK. Fowler was at 
this time already the leader of the world effort (which Caltech championed) to de- 
termine the rates at which nuclear reactions would occur in stars. This involved the 
now familiar technique of measuring nuclear interactions at MeV energies with Van 
de Graff accelerator beams of charged ions, and extrapolating measured data down- 



ward in energy to the Gamow-peak energy within a thermal distribution ( Clayton 



1968 Ch. 4). Fowler's zest for nuclear astrophysics was boundless, and he soon 
had the entire Kellogg Radiation Laboratory at work on thermonuclear reactions 
rates between positive ions. Fowler was awarded a share of the 1983 Nobel Prize in 
Physics for that pioneering program of research. His efforts had been directed pri- 
marily to the question of thermonuclear power in stars; but Fowler became intensely 
attracted also to the larger question of the origin of the elements. In Cambridge he 
met Geoffrey and Margaret Burbidge, and these three began a project with Fred 
Hoyle to write a survey paper on the issues of nucleosynthesis in stars. That pa- 



per was written at Caltech in 1956 and published the next year (Burbidge et al. 
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1957 i. Soon called simply B FH, it became one of the most celebrated papers in 



astrophysics. 

The authors of B 2 FH contributed creatively and energetically to formulating the 
neutron-capture processes for synthesizing the elements heavier than nickel. Ascrib- 
ing crucial roles to the magic neutron numbers N = 50, 82 and 126 they described 
environmental situations in stars within which selected isotopes would be abun- 
dant, fleshing out ideas that Suess had envisioned earlier for two neutron-capture 
processes. Suess had, however, been unable to formulate these clearly enough to 
win contemporary acclaim. The slow neutron-capture process, named s process by 
B 2 FH, envisioned 100-to-1000 years between neutron captures, so that radioactive 
isotopes would generally beta decay before capturing another neutron, keeping the 
capture path trailing along the valley of beta stability. The rapid-neutron-capture 
process, named r process by B 2 FH, envisioned neutron densities so large that neu- 
tron captures occur in tens of milliseconds, faster than beta decays, with the result 
that the capture path moved into the realm of radioactive neutron-rich nuclei, be- 
ing halted only when no additional neutrons could be stably added owing to their 
diminishing separation energies. These created waiting points at which the capture 



flow would halt and wait until beta decay occurs (see also Seeger et al. 1965 Clay- 



ton 



1968 ). Defining the s process and the r process was the high point of the B 2 FH 
contributions to nucleosynthesis theory. Their Appendix included an inspiring ta- 
ble of all heavy isotopes in which each was characterized as being either s-process 
only, mostly s process, comparable s process and r process, mostly r process and 
r-process only. This can be regarded as the next important step in nucleosynthesis 
theory. 

An important aspect of astronomy with radioactivity lies in the competition be- 
tween beta decay and neutron capture that ensues when neutron capture by stable 
isotopes create isotopes that can undergo beta decay. Most do so quickly in compar- 
ison with the time required to capture another neutron, but some key branch points 
are slow to decay. B 2 FH had inventively shown that such competition at branch 
points could, when compared with the actual solar abundances, reveal the time scale 
and neutron density for the operation of the s process. Thermally populated excited 
states of radioactive nuclei often increase the effective decay rate, a delicate point 
in that aspect of astronomy with radioactivity. 

Two cautions about B 2 FH must be made in order to not overly eulogize what 
they achieved. Firstly, both the s process and the r process are, as described by 
them, secondary processes of nucleosynthesis. Secondary nucleosynthesis refers to 
the synthesis of new heavy nuclei from other existing heavy nuclei. It does not in- 
crease the galactic metallicity, a goal that Hoyle's founding papers had achieved 
spectacularly as the primary goal of the astronomy of nucleosynthesis. The B 2 FH 
neutron-capture processes instead change one existing heavy nucleus into another. 
So those processes did not contribute to increasing metal abundances in the galaxy. 
Today it is known that, contrary to their description, the r process is actually pri- 
mary, because the collapsed supernova core synthesizes the seed nuclei that rapidly 
capture its free neutrons. The details of this are still not understood, however. As- 
tronomical observations of old metal-poor stars confirm that the r process indeed 
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began earlier than did the s process. Secondly, the B 2 FH descriptions of the s pro- 
cess and the r process could not immediately be used for astrophysics calculations 
because they were time-independent formulations. Both the s process and the r pro- 
cess were described by the static condition dN/dt = 0. The B 2 FH descriptions were 
thus of stationary abundances that could exist within appropriate environments for 
these processes. This enabled a rough but clear correlation between the nuclear sys- 
tematics of their process abundances and the solar isotopic abundances. But it did 
not enable calculation of the temporal growth of these abundances. 

Over the next two decades B 2 FH nonetheless became the default citation for 
workers wanting a reference to the general theory of nucleosynthesis in stars, vastly 
eclipsing the rate of citations to Hoyle's previous papers. One key to its success was 
its citation of over 100 stars showing nucleosynthesis effects in their spectra and of 
more than 100 astronomical research papers concerning those stars. This drew the 
astronomical community into the scientific culture of nucleosynthesis (which was 
then new to astronomers) as Hoyle's papers could not. Astronomers generally cited 
B 2 FH rather than Hoyle, with the result that Hoyle's papers slipped into relative 
obscurity. 

Fowler himself coauthored several other important works for nucleosynthesis, 
primarily with Fred Hoyle on supernovae and with Donald Clayton on nucleosyn- 
thesis by neutron capture chains; however, his most important and essential role in 
history was the empirical thermonuclear reaction rates determined in his laboratory. 
For this he shared the 1983 Nobel Prize in physics. Probably Fowler's most signif- 
icant subsequent work with Hoyle concerned the puzzling nature of the supernova 



phenomenon and on nucleosynthesis within them (Hoyle et al. 1960| l. They divided 
supernovae into Types I and II based on theory, rather than on the presence of H ab- 
sorption lines in the spectra. Type I was assigned to low-mass progenitor stars that 
evolved to white-dwarf stars of degenerate carbon, but then ultimately explode in 
an exothermic thermonuclear display that results in most of the iron in the universe. 
Electron-degenerate matter, supported by degeneracy pressure, is violently unstable 
to thermonuclear runaway. The Type II supernovae, on the other hand, occurs in 
massive stars whose cores are too massive to form white dwarf structure. In these 
continuing nuclear burning eventually exhausts the nuclear energy supply, and so 
the core must collapse. Fred Hoyle had predicted this inevitable collapse in 1954; 
but their 1960 paper carried the physical picture further. Today the world of astron- 
omy uses this Hoyle-Fowler classification based on the physics of the device. It was 



a landmark in astrophysics (see Woosley and Weaver 1999 1. 



A. G. W. Cameron, a Parallel Force 

The year 1957 also saw the emergence of A. G. W. Cameron as one of the pio- 
neers of nucleosynthesis. The Atomic Energy of Canada at Chalk River published 
in bound-mimeograph form a series of lectures on nucleosynthesis that Cameron de- 
livered at Purdue University in March/ April 1957 (Came ron] |1957| l. These lectures 
covered in an independent way essentially the same material as in B 2 FH. Cameron 
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had constructed his treatments by working alone, first as a new faculty member at 
Iowa State University and then at Chalk River, with only Hoyle's two papers to guide 
him as well as his training in nuclear physics from University of Saskatchewan. 
Owing to his independent treatment, Cameron's Lecture Notes became a valuable 
source of new ideas in nucleosynthesis. In his emphasis on nucleosynthesis within 
the separate shells in massive stars, Cameron's approach followed Hoyle's equa- 
tion and thereby enlarged the ideas in Hoyle et al. ( 1954) ) more effectively than did 
B 2 FH. For the next two decades Cameron stressed in many publications with re- 
search students - ones that he had recruited during a guest lecture course at Yale 
University - the vast nucleosynthesis changes that occur in a massive star at its time 
of explosion. His papers establish him as the first great disciple of Hoyle's equation. 
Cameron's later Yale Lecture Notes (1963), written by his Yale students W.D. Ar- 
nett, C.J. Hansen and J.W. Truran, were much improved over Cameron ( 1957). They 
probably should have been published, but were not and are therefore not generally 
available; however, they were a strong influence on history, especially on Cameron's 
students and on Clayton and his students. Partly owing to their unavailability, Clay- 
ton published his own textbook ( |Clay ton) 1 1 968) 1 . 

Cameron began his nucleosynthesis research, however, with the sources of free 
neutrons in stars that could be responsible for the observed radioactive Technetium 
observed by Paul Merrill in stellar atmospheres. Because all isotopes of the element 
Tc are radioactive, the presence of its absorption lines in stellar spectra argued that 
it had been created within the star during roughly its last lifetime. In 1955 Cameron 
proposed that the 13 C(a,n) 16 reaction would liberate the extra bound neutron in 
the 13 C nucleus and that the liberated neutron could be captured to create heavier 
isotopes ( |Cameron 1955) >. Later Cameron ( 1959) > calculated with the same moti- 
vation the number of neutrons liberated during carbon thermonuclear reactions fol- 
lowing the exhaustion of He. The carbon fusion reactions were part of the sequence 
of thermonuclear stages within evolving massive stars (as |Hoyle et al.| ([1954} had 
first described). 

In his large subsequent body of work, Cameron established himself as a true 
polymath. He introduced the speeding up of beta decay rate by thermal popula- 
tion of the excited states of a nucleus, many of which decay more rapidly than the 
ground state. During 1955-56 Cameron introduced numerical computation on the 
first vacuum-tube computers into nucleosynthesis problems, and remained there- 
after on the cutting edge of nucleosynthesis computation. Cameron and his students 
repeatedly blazed new paths by programming nuclear reaction networks onto the 
latest and newest computers. He became an expert on planetary sciences and an im- 
portant advisor to NASA. In particular, Cameron is forever famous for his work on 
the origin of the moon as the result of a Mars-like planetary collision with the young 
earth, a theory of origin that now seems beyond doubt. 
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Donald D. Clayton and Time-Dependent Heavy Element Nucleosynthesis 

Beginning as research student at Caltech with Fowler, Donald Clayton began con- 
structing a time-dependent formulation of the s process in 1957, the year of B 2 FH 
publication and of Cameron's Chalk River lecture notes. Clayton's discard of the 
assumption dN/dt = 0, an assumption requiring a constant aN curve for s-process 
abundances, altered profoundly the direction of s-process research by focusing on 
how efficiently seed nuclei could be converted to heavy s-process nuclei. B 2 FH had 
not addressed that question. Clayton showed, as B 2 FH had surmised, that the iron 
abundance peak must provide the seed nuclei being transmuted into the large over- 
abundances of barium in stars whose spectra showed Ba/Fe some 20-50 times the 
solar ratio. More surprisingly his results also showed that, as cumulative neutron 
fiuence increases, none of the sequential abundance distributions that are generated 



resemble the solar abundances ( |Clayton et al. 1961 1. The solar s-process abun 



dances were required to be a superposition of differing numbers of Fe-seed nuclei 
(per Fe nucleus) exposed to differing integrated fluxes of free neutrons. The number 
exposed must be increasingly smaller for increasingly larger neutron irradiations. 
Galactic history or history within s-process stars is required to bring that superpo- 
sition about. Therein lay new astrophysics. The solar s-process abundances were 
shown to not be simply a smoothly declining aN curve, as B 2 FH had speculated, 
but a superposition of exposures generating narrow regions of atomic weight near 
the neutron-magic numbers where the assumption dN/dt = is severely violated. 
Two decades of improved measurements and consequent fitting to solar abundances 



( Kappeler et al. 1982| l would be required before advances in stellar evolution would 



be able to describe the exponential-like fiuence distribution that was required. It was 
a sophisticated interplay between He-shell pulses and cyclically ignited H burning 
at the base of the envelope of AGB stars (see Chapter 3). 

From the time of these first solutions of the neutron-irradiation superpositions 
resulting in the s-process abundances, new phenomenological aspects of heavy el- 
ement nucleosynthesis were possible. The theory-based fit yielded all s-process 
abundances with meaningful accuracy. These allowed Clayton and Fowler (1961) 



to publish the first decomposition of heavy-element abundances into their s-process 
and r-process parts. B 2 FH had suggested the dominant processes for each isotope in 
their spectacular appendix; but a quantitative decomposition became possible for the 
first time. Their initial effort has been redone at least a dozen times as new neutron- 
capture-cross-section data appeared, most notably first by See ger et al.| ( |1965] ) and 
later by a new measurement program in Karlsruhe (Kappele r et al.| 1 1982 ). This 
s-r decomposition applied to astronomical spectroscopy of old stars has routinely 
produced meaningful new knowledge. Observations of old metal-poor dwarf stars 
indicated that the r-process abundances began to grow earlier as star formation first 
began than did the s-process abundances (Truranj 1981 Burris et al. 2000). That 



result demonstrated that the r process is a primary nucleosynthesis process, rather 
than secondary as B 2 FH had stated. This requires the r-process to occur within core- 
collapse supernovae. Quantitative s-r decomposition also inspired unforeseen new 
techniques for radioactivity-based cosmochronology. 
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Impressed by the new astrophysics lurking in time-dependence, Clayton advo- 
cated time-dependent formulation of the r process as well, again jointly with Fowler 
and with P. A. Seeger, Fowler's research student in Kellogg Lab. For the r process 
as for the s process, B 2 FH had described only a time-independent steady flow that 
showed the neutron-rich heavy isotopes did indeed have abundant progenitors in 
suitable time-independent settings; but they had not been able to address whether 
the entire r-process abundances can be synthesized at one set of conditions. Com- 
puters at that time (1963) were not capable of handling a full r-process network. 
Seeger et al. (1965 ) showed that the full mass range cannot be produced together 



(unless new seed nuclei are injected during the process). B 2 FH had creatively de- 
fined the key nuclear physics relationships of the r process but were mute on its 
dynamics. The time-dependent formulation by Seeger et al7|( |1965| > became a proto- 
type for r process astrophysics. It showed the r process to also be a superposition of 
differing irradiation histories of seed nuclei. Even after four decades of subsequent 
computations, the nature of the r process superposition remains a frontier puzzle. 
It is noteworthy for astronomy with radioactivity that the entire r process reactions 
occur within the realm of radioactive nuclei. Only after rapid expansive cooling can 
that neutron-rich radioactive abundance distribution undergo a series of beta decays 
changing each isobar identity until resting at the most-neutron-rich stable isobar 
(isobar is an isotope having the same atomic weight). A major research goal of nu- 
clear astrophysics today is better laboratory definition of the parameters defining 
accurately the properties of the neutron-rich radioactive nuclei. 

Understanding the time dependence during silicon burning (silicon photoero- 
sion) was the big scientific challenge of the mid-1960s. Clayton introduced nuclear 
quasiequilibrium as a physical concept (Bodansky et al. 1968a|b| l to clarify how sil- 
icon transmutes to an iron abundance peak. Quasiequilibrium explained the only big 
gap in Hoyle's 1954 theory of primary nucleosynthesis, replacing the ill-formulated 
a process of B 2 FH. A temporal sequence of quasiequilibrium states facilitated the 
calculation of the set Am^ for k=28-62 to be inserted in Hoyle's equation. The se- 
quence of quasiequilibria again involved relaxing the assumption dN/dt = 0. The 
quasiequilibrium concept was powerful and new, and enriched many subsequent 
aspects of nucleosynthesis reaction networks (the r process, explosive oxygen burn- 
ing, the a -rich freeze out, the origin of 48 Ca, and others). Just as all nuclear reac- 
tions proceed at the same rate as their inverses in full nuclear equilibrium (NSE), 
during quasiequilibrium one refractory nucleus violates that equilibrium assump- 
tion by changing abundance only slowly, while all others maintain equilibrium with 
it. During silicon burning the 28 Si nucleus is the slowly changing, refractory post 
to which the quasiequilibrium distribution is attached (B odansky et al.| [l968a|b[ l. 
The most abundant isotopes between A=44 and 62 are, in this quasiequilibrium 
sequence, created as radioactive progenitors rather than as stable isobars, with im- 
portant consequences for emerging astronomies of radioactivity. 
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The Sequel 

The decade 1956-66 following Hoyle's pioneering two papers had witnessed pro- 
found enlargements of the theory of nucleosynthesis in stars. The years following 
the publication of B 2 FH had been marked by vast improvement and reformulation 
of its influential processes. Calculable time-dependent descriptions of heavt-element 
nucleosynthesis processes refocused attention from simple correlations between nu- 
clear properties and abundances to the astrophysical histories and stellar evolution 
that bring them into existence. The sometimes heard statement by astronomers that 
not much happened after B 2 FH reflects lack of awareness of these historic changes. 
It may be of interest to note that the beginnings of nucleosynthesis theory was an in- 
ternational innovation. Of the pioneers named, only Fowler and Clayton were born 
Americans; Hoyle, E.M. Burbidge and G. R. Burbidge were English; Suess was 
Austrian; and Cameron was Canadian. 

Innovations continued and accelerated during the next decade 1967-77. These 
will not be reviewed here except to say that the evaluation of Hoyle's equation 
through numerical computation of the evolution of massive stars yielded repeated 
insights into the interplay between stellar evolution and nucleosynthesis. The B 2 FH 
neutron-capture processes revealed layers of complexity associated with the stars 
and with time dependences. The innovative center of this research moved away 
from Caltech, initially to Yale University and to Rice University, where Cameron 
and Clayton respectively founded schools evaluating Hoyle's equation. New leaders 
developed within those schools included especially W. David Arnett from Yale (sub- 
sequently also Rice) and Stanford E. Woosley from Rice (subsequently U.C. Santa 
Cruz). Each has many important publications on the evolution of massive stars and 
the nucleosynthesis in its shells ( |Arnett 1996 [Wbosley and Weaver 1995) (see 



Chapter 4). A European role in this research also experienced rebirth during that 
decade, especially in Munich. It is also the decade 1967-77 that sees the emergence 
of several new observational aspects of astronomy with radioactivity. A description 
of those developments follows. 



2.1.3 New Astronomy with Radioactivity 

The existence of natural radioactivity clearly holds implications for the origins of 
atomic nuclei. The chemical elements could not have always existed if the radioac- 
tive nuclei were created along with the stable nuclei. In this way naturally occurring 
radioactivity is intimately related to nucleosynthesis of the chemical elements. It 
was evident from the structure of the suggested processes of nucleosynthesis that 
radioactive nuclei played a large role in each of them. The key role of radioac- 
tive isotopes during nucleosynthesis and during thermonuclear power in stars was 
the first astronomy with radioactivity. Presumably the radioactive nuclei would be 
ejected from stars along with the new stable nuclei unless they decayed within the 
stars prior to ejection. Radioactivity plays a large role in the H-burning reactions, 
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the PP chains and the CNO cycle (see |Clayton| |1968| l responsible for the stellar 
power capable of keeping the stars from cooling. Hans Bethe was the 1968 Nobel 
Prize awardee for discovering these H-burning cycles in stars just prior to WWII. 
So it was evident in mid-20 th centrury that radioactive nuclei carry significant issues 
for astronomy. In the 1950s the idea arose of viewing the radioactivity at the solar 
center by detecting neutrinos arriving from the sun. Raymond Davis Jr. would win 
a Nobel Prize for spearheading that effort. 

What made astronomy with radioactivity so exciting scientifically was the dis- 
covery in the 1960s and 1970s of altogether new ways of observing radioactivity in 
astronomy. These aspects of astronomy with radioactivity generated new interdis- 
ciplinary connections to nucleosynthesis. All science needs observations to provide 
an empirical base, astronomy just as surely as laboratory science. To be sure, the 
scientific method happens differently in astronomy than in laboratory science. In 
astronomy there are no experiments that can change the initial conditions as a test 
of theories. No experiments can be planned to refute a hypothesis. One has instead 
only observations of natural events. In many cases simple observations can be suf- 
ficient to refute a hypothesis. It is fortunate that nature provides so many natural 
events, so that in many cases contrasting separate events constitutes a type of ex- 
perimentation, giving a spectrum of observations in which conditions differ in ways 
that must be inferred but are nonetheless real. The experiments are natural and per- 
formed by nature herself, however, rather than by scientists. 

New types of observations of the occurrence of natural radioactivity galvanized 
astronomy by providing new kinds of astronomical data. These first historical ad- 
vances are reviewed in what follows, as previews to the more detailed chapters that 
discuss them. They are: cosmochronology, the age of the elements; gamma-ray-line 
astronomy of radioactivity; radioactive power for light curves of supernovae and of 
novae; extinct radioactivity measured by excess abundances of daughter nuclei in 
solids formed in the early solar system; extinct radioactivity measured by excess 
abundances of daughter nuclei in solids condensed while dust particles condense as 
gas leaves a site of stellar nucleosynthesis. What follows is a brief description of 
how these applications were discovered or anticipated. 



2.1.4 Interpreting Exponential Decay 

The familiar exponential law of radioactive decay follows from Eq. (1.1) for an 
ensemble of N radioactive nuclei from the constancy of decay probablility for a 
single nucleus. Given N such nuclei the expected number of decays per unit time 
in the ensemble is given by the product of the number N of nuclei and the decay 
probablility per unit time for a single nucleus. Accordingly 



(2.4) 
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Because A, and therefore z, is a constant for that nuclear species when it exists 
outside of stars, this equation has a well-known integral form, 

N(t)=N e-'/ T (2.5) 

where No is the initial number in the ensemble (at f=0). Similarly, if a number A^o is 
observed now at t = to, the number that existed at an earlier time t \ would have been 

N(t l )=Noe i '°-'< )/T (2.6) 

provided that new nuclei have not been added to the ensemble during that time 
interval. In the same spirit, if N\ were created at earlier time t \ and N2 were added 
at a subsequent time t2, the number A^o that would exist today at to is 

^0=^1 e-^-' l)lx +N 2 e- (, °- ,2) l % . (2.7) 

Generalizing slightly, let B & (t) represent the fractional age distribution of the 
primary stable solar nuclei at the time of solar formation at t — t e . Then dB & {t)/dt 
is the number of stable solar nuclei that had been born per unit time at time t . It is 
the age distribution that a radioactive species within an ensemble would have if it 
were not decaying. The age distribution for radioactive parents may be thought of as 
the age distribution of stable nuclei that were synthesized along with the radioactive 
nuclei. A mnemonic for the symbol Z? fl is the birthrate of those solar system nuclei. 
dB, 3 jdt is the rate at which the stable nuclei and the primary parent radioactive nu- 
clei were added to the total sample destined for the solar system rather than the rate 
at which they were produced by nucleosynthesis. In drawing this subtle distinction, 



Clayton ( 1988) demonstrated the error of the common practice of equating the age 



distribution of solar nuclei with the rate of galactic nucleosynthesis. The age dis- 
tribution of solar-system nuclei is vastly different from the galactic nucleosynthesis 
rate because the evolution of the ISM is involved. Even if the ISM is instantaneously 
mixed, many of the old nuclei became trapped within stars, so that the age distribu- 
tion in the ISM becomes biased toward more recent nucleosynthesis. 

Given this definition of dB & /dt as the age distribution of solar-system nuclei, 
and given that their addition to the solar sample began at time t\ and continued until 
?0, the total number of any stable primary nucleosynthesis product would have been, 

f'& dB G) (t') , 
N tot ai = J t -^y- dt' (2.8) 

Then in analogy with Eq. |2.7| if those nuclei are instead radioactive, the number 
surviving until solar formation at time f would be 



[ t@ ^ e-^-OIW (2.9) 
Jti dt' 

If the total duration to — f 1 of production is much greater than the mean lifetime T, 



and if dB^/dt varies only slowly, Eq. 2.9 reduces approximately to 
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This useful formula estimates the number of remaining radioactive nuclei to be equal 
to the number created during the last mean lifetime of that radioactive species. If 
dB {?} /dt has jagged temporal structure near f , however, a more complicated evalu- 
ation would be required. The Eq. |2.10| is very useful as a first approximation to the 
numbers of radioactive nuclei within the ISM during continuous nucleosynthesis in 
stars if the birthrate B{t) is redefined as the birthrate spectrum of those stable ISM 
nuclei residing in the ISM rather than in stars. 

These relationships expressing properties of exponential decay are central to un- 
derstanding both extinct radioactivity in the early solar system and the numbers of 
radioactive nuclei within astronomical objects. 



2.2 Disciplines of Astronomy with Radioactivity 
2.2.1 Nuclear Cosmochronology 

Today it is self evident that the existence of radioactive nuclei implies that they 
were created at some estimable moment in the past. Were that not the case, they 
would have long ago have disappeared. But the full implications were not evident 
to those engaged in the exciting day-to-day goals of understanding natural radioac- 
tivity. Nonetheless, in 1929 Ernest Rutherford wrote what may be the first paper 
on astronomy with radioactivity (Rutherford 1929| >. Today we use the term nuclear 



cosmochronology to mean the attempt to use natural abundances of radioactive nu- 
clei to compute the age of an astronomical object or of the elements themselves. 
Determining the age of the chemical elements by assuming them to have the same 
age as the radioactive nuclei became a goal that has attracted many, most notably 
William A. Fowler and this writer. 



Uranium and Thorium on Earth 

When he began to think on these things, Ernest Rutherford concluded from the 
ratios of 235 u/ 238 U as measured in his laboratory that uranium was created somehow 
within the sun and transported to earth. Accepting Jeans's estimate that the age of 
the sun was an exuberant 7x 10 12 years, he wrote in a prophetic paper addressing 



the puzzle (Rutherford 1929): 



..it is clear that the uranium isotopes which we observe on earth must have been forming in 
the sun at a late period of its history, namely, about 4xl0 9 years ago. If the uranium could 
only be formed under special conditions in the early history of the sun, the actino-uranium 
on account of its average shorter life would have practically disappeared long ago. We may 
thus conclude, I think with some confidence, that the processes of production of elements 
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like uranium were certainly taking place in the sun 4xlCP years ago and probably still 
continue today. 

Corrected modern numbers for those that Rutherford used are the observed 
abundance ratio 235 U/ 238 U=7.25x 10~ 3 and the mean lifetimes against alpha de- 
cay, T( 235 U)= 1.029xl0 9 yr and t( 238 U)= 6.506x10 9 yr. Since each isotopic 
abundance has been exponentially decaying during the age An of the earth, their 
abundance ratio when earth formed would have been ( 235 U/ 238 U)o = 7.25 xlO~ 3 

(e A E /T( 235 U) /^ £ /T( 235 U)^ 

The age of the earth has been reliably measured using the fact that these two 
isotopes of U come to rest, after a series of alpha decays, as different isotopes of Pb, 
namely 207 Pb and 206 Pb respectively. Using the measured ear thage, A E = 4.57 xlO 9 
yr, yields the initial U ratio on earth to have been ( 235 U/ 238 U)o = 0.31. These facts 
are beyond doubt. 

The implication of great consequence stems from the expectation that the r pro- 
cess that is responsible for the nucleosynthesis of both isotopes, should make more 
235 U than 238 U. 235 U has six non-fissioning progenitors whereas 238 U has but three. 
An r-process progenitor is a non-fissioning transuranic nucleus that after a series 
of relatively fast radioactive decays comes to rest at one of these long-lived iso- 
topes of U. For 235 U those nuclei are A =235, 239, 243, 247, 251 and 255-totaling 
six r-process progenitors-whereas for 238 U they are A= 238, 242, 246 and 35% 
of 250-totaling 3.35 progenitors. Taking into account the empirical evidence that 
production by the r-process favors even-A nuclei by a slight 20% over adjacent 
odd-A nuclei, one expects 235 U production to exceed that of 238 U by a factor near 
P( 235 U)/f( 238 U)=1.79. This argument reveals that the abundance ratio 235 U/ 238 U 
has declined from near 1.79 at production to 0.31 when the earth formed. This de- 
cline takes considerable time, showing that U isotopes were synthesized during pre- 
earth astrophysical history. 

This first calculation of nuclear cosmochronology reveals the nature of the prob- 
lem, but also its uncertainties. Is the production ratio f( 235 U)/P( 238 U)=1.79 cor- 
rectly estimated? Were the U isotopes synthesized in one single presolar event, in 
which case it occurred about 6.6 billion years ago, 2 billion years prior to formation 
of the earth. Or was their production rate distributed in presolar time? If the age 
distribution of the solar system's r-process nuclei is flat between the time of first 
production and the earth's formation, the beginning of nucleosynthesis would fall 
near 13 billion years ago. The true age distribution can hardly be known with any 
assurance, revealing the severe limitation of this single pair for constraining the time 
of the beginning of nucleosynthesis in our galaxy. 



The technique was extended to the ratio 232 Th/ 238 U by Fowler et al. (1960); 
and Fowler returned to it many times in later years (e.g. [Fowler et al. 1972| l; see 



also the textbook by Clayton (T968J). Relevant numbers used are the observed abun- 



dance ratio 232 Th/ 238 U= 4.0 and the mean lifetimes against alpha decay, r( 232 Th)= 
20.04x 10 9 yr and t( 238 U)= 6.506x 10 9 yr. The abundance ratio is much more un- 
certain than it was for the pair 235 U/ 238 U because Th and U are different chemical 
elements. Any two elements having different fractionation chemistry pose a tough 
problem when seeking their ratio in the initial solar system (or in the sun). Since 
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each isotopic abundance has been exponentially decaying during the age Ae of the 
earth, their abundance ratio when earth formed would have been ( 232 Th/ 238 U)0=2.5 
instead of 4.0. And the production ratio in r-process events was inferred from the 
same counting arguments of odd and even progenitors of 232 Th and 238 U to be 
P( 232 Th)//>( 238 U) = 1.73. 

A problem is that this Th/U pair does not give transparently concordant num- 
bers with the U pair. Uncertainties in r-process production ratios, in the relative 
abundances of Th and U, and in the arbitrary parameterizations of galactic chemi- 
cal evolution that have been used each conspire to yield possible solutions in which 
nucleosynthesis began anywhere from 2 to 10 Gyr prior to solar birth. Fowler strove 
repeatedly to circumvent these uncertainties, trying to extract the true answer (for 
there obviously is a true answer!); but in fact, the data are not adequate and the as- 
trophysical model of the galaxy used is not adequate for the task. Data from other 
chronological species would be needed, along with a more sophisticated apprecia- 
tion of galactic chemical abundance evolution. 

On a more positive note, Fowler's papers (Fowler etal. 1960|b 1972 and several 
others) inspired many others to tackle this fascinating aspect of astronomy with 
radioactivity. 



Cosmochronology with Extinct Radioactive Abundances 



In 1960s a new discovery, excess trapped 129 Xe gas in meteorites (Reynolds et al 



1960), provided the first evidence that the matter from which the solar system 
formed contained radioactive nuclei whose halflives are too short to be able to sur- 
vive from that time until today. These are called extinct radioactivity. Discovery of 
extinct radioactivity utilized the buildup of the daughter abundance from a radioac- 
tive decay as a mea sure of how abundant t hat radioactive parent was in the initial 
solar system. What Reynolds et al. ( 19601 showed was that the excess 129 Xe gas 
in meteorites had resulted from 129 I decay and that the initial 129 I abundance was 
about 10~ 4 of the initial 127 I abundance. He realized that this datum allowed an es- 
timate of when the iodine isotopes had been created. By assuming that at the time 
of their nucleosynthesis the initial 129 I abundance had been equal to that of stable 
initial 127 I, Reynolds argued from the mean lifetime t( 129 I)= 23.5x10 6 yr (23.5 
My) that iodine had been created only about 300 My prior to solar system birth. 
This was the first new cosmochronology technique since Rutherford and brought 
Reynolds great fame. Unfortunately, Reynolds' conclusion was no more believable 
than Rutherford's had been, because his astrophysical model was also quite unre- 
alistic. He mistakenly assumed that all of the iodine abundance in the solar system 



had been created at the same time, and therefore used Eq. 2.5 to estimate that time. 
Because 10~ 4 requires a decay period equal to about 13 mean lifetimes, he had con- 
cluded that the age of the elements was 13t( 129 I) 



Fowler et al. 



(1960b) 



300 My. 

countered that it was more realistic to assume that the synthesis of iodine isotopes 
was spread out by multiple r-process events, occurring between the time ti of a first 
nucleosynthesis event until the time ?o of the last such event prior to solar birth. That 
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concept, called continuous nucleosynthesis, required use of the Eq. 2.10 limit of Eq. 



2.9 instead of Eq. 2.5 That calculation gave the initial nucleosynthesis epoch to be 
104t( 129 I) = 230 Gy - very old indeed, unrealistically old. This picture was then 
modified to allow for an interval of sequestered interstellar gas, probably within a 
molecular cloud, of about 100 My during which no new nuclei were added to the 
solar mix - a so-called free decay interval in which radioactive abundances would 
follow Eq. 2.5 Such an interval would require that 100 My before the solar birth, 
when that interstellar gas withdrew from new nucleosynthesis, the interstellar abun 
dance ratio would have been /V( 129 I)/Ar( 127 I) = 2.8 x 10~ 3 rather than 10~ 4 . Eq 



2.10 



then suggests that nucleosynthesis began 8 Gy prior to solar birth, a reasonable num- 
ber, obtained, unfortunately, by construction. But that construction had nonetheless 
established a new paradigm for radioactive abundances in the early solar system. 
The upshot for the problem of radioactive chronology is that the abundance of now 
extinct radioactivity in the early solar system can say little about when nucleosyn- 
thesis first began. But it provides other equally interesting issues for the astronomy 
of radioactivity. 

Radioactive 129 I is of special historical significance as the first extinct radioactive 
nucleus to be discovered that can be regarded as part of the average radioactivity in 
the galactic ISM at the time of solar birth. If so it is in that sense typical of what 
is expected in the ISM. Other extinct radioactivities may not be capable of such 
an interpretation but must instead be interpreted as produced by a special event 
associated with the birth of the solar system. Other subsequently discovered extinct 
nuclei that can be regarded in the first category of galactic survivors are 
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Hf, 146 Sm and 244 Pu. 



Birthrate Function for Primary Solar Abundances 

The history of nuclear cosmochronology reveals widespread confusion between the 
rate of galactic nucleosynthesis and the age spectrum of solar system nuclei. It is 



the latter that enters into cosmochronology through Eq. 2.9 So endemic is this con- 
fusion between these two fundamental concepts that it is important to define in this 
section the birthrate spectrum for primary abundances in the solar system. 

Let Nq be the solar abundance of a primary nucleosynthesis product. One may 
think of it either as the total number of those atoms in the solar system, which 
consists of awkwardly large numbers, or as the total number of those atoms in the 
solar system normalized to a defined abundance for a specified nucleus that sets the 
abundance scale. Two such normalized scales are very common: one used more in 
astronomy defines the solar abundance of hydrogen to be Nh = 10 12 ; another, used 
more in solar abundances derived from meteorites, defines the solar abundance of 
silicon as A^,=10 6 . This book will primarily use the latter definition because most 
studies of nucleosynthesis in the literature use the scale Si=10 6 (where the chemical 
symbol is often used to represent the abundance of that element). Let B Q (t) be 
the cumulative fraction, as a function of time, of stable solar primary nuclei that 
already existed at time t. Letting for convenience the starting time t = be the time 
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when nucleosynthesis first began, the function .60(f) starts at at f=0 and rises to 
Bq=\ at t = to, the time when the solar system formed. Note carefully, B & (t) is not 
the cumulative birth fraction of that nucleus in galactic nucleosynthesis, but rather 
the cumulative birth fraction of those primary nuclei that actually entered the solar 
system inventory. Normalizing this definition of #©(/) as a function that rises to 
unity at solar birth, the product 



N e -B @ (t)=N(t) (2.11) 

where Nq is the solar abundance of that stable species, N(t ) is the number of of those 
solar-system nuclei that already existed at time t . It is equally clear that dB & (t)/dt is 
the fractional birthrate per unit time of primary solar nuclei. It is the age distribution 
of solar nuclei, and redefined as here to rise to unity it is approximately the same 
function for all primary stable nuclei. 

A mental experiment is needed to set this birthrate function clearly in mind before 
considering why it differs so from the nucleosynthesis rate. The reader is warned 
that the research literature suffers endemic confusion over that distinction, requiring 
its careful definition. For any primary nucleus within solar abundances, paint each 
solar-system atom of it red! This cannot of course be done in reality, but the thought 
experiment helps understanding of the birthrate spectrum. On the scale Si=10 6 one 
has by construction 10 6 red Si atoms in the solar system. Then imagine watching 
a film run backwards in time of these 10 6 red Si atoms. Back in the interstellar 
medium they are mixed with a much larger number of unpainted Si atoms. As the 
film runs to earlier times, these 10 6 red Si atoms have greatly differing histories. 
Some, after being created in a supernova, have later been inside a star, survived and 
come back out. But for each atom, its world line reaches a time when that atom first 
came into existence. Each of the 10 6 red solar Si atoms is labeled by that birthdate, 
born from distinct supernovae at differing times. Then a graph is constructed giving 
the total number that existed at time t . This graph starts at 1 atom when the first solar 
Si atom is created and reaches 10 6 by the time the solar system forms, when all 10 6 
have been created. If the curve is then normalized by dividing by 10 6 , one has the 
birthrate spectrum B & (t) of solar primary nuclei, which rises to unity at the time 



the solar system forms. Eq. 2.11 depicts this situation, where A^(f) rises to 10 6 at 
solar birth, but B, 3 (t) rises to unity at that time. Although different primary species 
have vastly differing abundances, for each the function B^(t) is approximately the 
same function because the abundances of all primary nucleosynthesis products rise 
approximately (but not exactly) together. 



Eq. 2.11 is the solution of the differential equation 

dN/dt = dB & (t) jdt for stable nuclei (2. 12) 



and 



dN/dt = dB, ?) (t)/dt — XN for radioactive nuclei 



(2.13) 
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One frequently sees these equations written for galactic nucleosynthesis, with dB/dt 
identified as the galactic nucleosynthesis rate and N(t) as the ISM abundance. This 
is a fundamental error. Consider why. 

The galactic nucleosynthesis rate for Si atoms was very large in the young galaxy, 
when the rate of formation of massive stars and supernovae was large. But those new 
Si atoms in the interstellar medium become mostly locked up in the interiors of sub- 
sequent low-mass stars, and are not available for incorporation into the solar system. 
This occurs because it is the fate of most interstellar gas to be locked up in low-mass 
stars. The initially gaseous galaxy is today about 10% gas, with 90% trapped in 
low-mass stars. The origin of the Si atoms that entered solar matter is therefore bi- 
ased towards supernovae that occur relatively shortly before solar birth, and biased 
against those that were created long earlier. This important conceptual distinction 
was demonstrated quantitatively by Clayton (1988), who stressed the concept of 
the age spectrum of solar nuclei and who constructed analytic models of galactic 
chemical evolution that enable a full disclosure of such issues. As an example, the 
simplest of all realistic models supports a star formation rate and supernova rate 
that declines exponentially as the interstellar gas is exponentially consumed; but the 
age spectrum of solar nuclei in that model is constant - equal numbers from equal 
times, despite the strong bias of nucleosynthesis of Si toward early galactic times. 
The student can study Clayton ( 1988 1 for many related issues for cosmochronology. 
Some of these issues involving galactic chemical evolution will be addressed later. 



Uranium and Thorium in old Dwarf Stars 

An altogether new technique in radioactive chronology became possible as CCD 
detectors enabled astronomers to measure line strengths for much weaker lines than 
had been previously possible. But cher et al.| ( |1987] l advanced the first such argument 
when he was able to accurately measure the abundance of thorium in old dwarf stars. 
The measured abundances enabled him to argue that the oldest stars are no older than 
10 Gy. If they were, Butcher argued, the Th, which has resided in the old stars since 
their births, would by now have decayed to a smaller abundance than it is observed 
to have. Because the dwarf stars observed appear to be among the early stars formed 
in our galaxy, the argument concludes that the galaxy age is negligibly greater than 
10 Gy. It will be clear that since the 232 Th mean life is a very long 20 Gy, it can 
have decayed from its initial abundance by only e~° 5 = 0.6 during a 10 Gy life 
of a dwarf star. Therefore the method requires not only an accurate measurement 
of the Th abundance but also an accurate estimate of its initial abundance when 
the star formed. It is the reasoning required to obtain that initial abundance that 
is a controversial aspect of this method. Realize also that the Th in an old star is 
decaying freely, so that its abundance declines faster than in the ISM where new 
nucleosynthesis of Th continuously replenishes it. 

The method advocated measures the ratio of Th/Nd in the star and observes that 
the measured ratio is smaller than that seen in younger stars. It assumes that the 
initial Th/Nd ratio in the star would be the same as the ratio seen in young stars, 
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since both Th and Nd are products of the r process and may be hoped to have a 



constant production ratio there. It was then counter argued (Clayton 1987 1988) 



that a constant production ratio for Th/Nd is not to be expected because almost half 
of the solar Nd abundance has been created in the s process, and, furthermore, the 
r process is primary whereas the s process is secondary. There exists ample and 
exciting evidence from the observed r-process pattern of abundances in extremely 
metal-deficient old stars ( |Truran| [1981] |Gilroy et ai.||1988||Beers et all |20Q5| > that 
the r process nucleosynthesis began prior to the beginning of the s process. That 
seemed to support the initial skepticism about Butcher's technique. But astronomers 
hurried to point out that if one omits the very metal-deficient old stars, the ratio of 
r-to-i abundances in stars having more than 10% of solar abundances is observed 
to be a near constant. |Clayton| ( ] 1988[ ) revisited that larger puzzle of parallel growth 
of s and r nuclei, noting that the abundance evidence suggests that the s process is 
primary despite its building upon iron seed nuclei; moreover, he presented a nucle- 
osynthesis argument that showed how the s process could in fact resemble primary 
nucleosynthesis even though it is secondary! See section 3.1 of Clayton ( 1988 1 for 
that argument, which has proven to be of high significance for nucleosynthesis the- 
ory. 

The lesson to be taken from this exciting new technique and its controversy is 
this: radioactive nuclear cosmochronology is vitally dependent upon a correct pic- 
ture of the chemical evolution of the galaxy. As such it has become less the province 
of nuclear physics and more the province of astronomy. Tinsley and Clayton both 
made that point in timely and influential ways. Only when many essential details of 
the history of our galaxy and of the history of nucleosynthesis within it have been 
settled can these techniques of radioactive chronology yield a reliable answer for 
the age of the chemical elements. Nuclear cosmochronology is truly an astronomy 
with radioactivity. 



Cosmoradiogenic Chronologies 

In the early 1960s a different approach to radioactive chronology became possible. It 
utilizes the buildup of the daughter abundances of radioactive decay during the his- 
tory of interstellar matter as a measure of how long that decay had been occurring 
and, therefore, how long ago the production of radioactivity began. One can imagine 
the stable daughter of a radioactive decay as a bucket into which all interstellar de- 
cay of the radioactive parent has been collected. It integrates past decay rather than 
focusing on how much radioactivity remains. Such daughter-isotope buildup during 
the history of the earth was already known as radiogenic abundance. It had been 
applied to the ages of earth and of meteorites, samples in which it could be expected 
that the other initial isotopic compositions were well known-namely, the solar abun- 
dances. But radiogenic abundance collected during the history of interstellar matter, 
what Clayton ( 1964| l called cosmoradiogenic abundance, was not seen as possi- 



ble data for determining the age of the elements themselves because the interstellar 
abundances are increased by the processes of nucleosynthesis as well as by any ra- 
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dioactive decay for specific isotopes. Thus the cosmoradiogenic abundance could 
not be easily disentangled from the direct nucleosynthesis abundance. 

What cut through that impasse was a credible quantitative theory of heavy- 
element nucleosynthesis. Suess, B 2 FH and Cameron had each contributed impor- 
tantly to the idea that two distinct neutron-capture processes, the s process and the 
r process, had been responsible for the creation of all but the very lightest isotopes 
of the elements heavier than about Z=32. The disentanglement required two things: 
firstly, an accurate s-process theory and secondly, a parent-radioactivity abundance 
that could be produced only by the r process and a daughter that could be produced 
only by the s process. The radiogenic daughter of such a radioactive r-process iso- 
tope is called a shielded isotope, because it is shielded from r-process production 
by the radioactive parent. 

Clayton (1964) made the key first step by noticing that the solar abundance of 
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Os is about twice as great as it is expected to be from s-process nucleosynthesis, 
and that it cannot be synthesized by the r process because neutron-rich matter at 
A=187 will, upon decaying toward stability, arrest at 187 Re. The decay cannot reach 
187 Os, which is shielded by 187 Re. The r process production flows into 187 Re. The 
abundance of 187 Os can then be thought of as having two parts; a part produced 
by the s process and a part produced by the decay of very-long-lived 187 Re. The 



quantitative s-process theory (Clayton et al. 1961 1 reliably accounted for about half 
of the 187 Os abundance, so the other half had to be the result of the beta decay of 
187 Re during ISM residence: 

187 Re — > 187 Os + e~+v f 1/2 =43Gyr 

That the halflife exceeds the galaxy's age is useful, because only a modest fraction of 
the 187 Re can therefore have decayed during presolar history and most of it therefore 
still exists; but because 187 Re is six times more abundant than the s-process amount 
of 187 Os, a decay of only l/6 th of the 187 Re can have doubled the s-process 187 Os 
abundance. The persistence of live 187 Re is even necessary for the method, because 
it allows one to equate the total amount of 187 Re nucleosynthesis that was destined 
for solar incorporation to the sum of the quantity of 187 Re remaining at solar birth 
and the quantity of 187 Re that had already decayed prior to solar birth. The latter is 
the difference between the total solar abundance of daughter 187 Os and the quantity 
produced by the s process. That is, 

187 0s = 187 Os s + 187 Os c 

where 187 Os c is the cosmoradiogenic part of the 187 Os abundance owing to cos- 
mic 187 Re decay. And the ^-process part can be obtained from 187 Os s = 186 Os 
(c(186)/c7(187)), where a is the appropriate neutron-capture cross section during 
the operation of the s process. This last relationship is valid because the entirety of 
186 Os is from the s process because it is shielded from r-process nucleosynthesis. 

The discovery of the Re-Os clock was the key that Clayton] ( |1964| l used in pre- 
senting solutions to three cosmoradiogenic chronologies. These are the beta decay 
of 187 Re to 187 Os described above, the similar beta decay of 87 Rb to 87 Sr, and the 
a -decay chains by which 235 ^ 238 u decay to 2 , 206 Pb. Each had been known in the 
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study of meteorites, but by this work they joined studies of the age of the elements 
themselves. 

Cosmoradiogenic lead is the more interesting of the other two chronologies, 
in part because it couples cosmoradiogenic Pb to the older chronology based on 
235 U/ 238 U, and in part because the Pb isotopes have three contributions to their 
abundances: 

206p b _ 206p bj + 206p br + 206p bc 
207p b = 207p bv + 207p br + 207p bt 

where 206 Pb c is the cosmoradiogenic part from 238 U decay, and where 207 Pb c is the 
cosmoradiogenic part from 235 U decay. Suffice it to be said here that the s-process 
part and the r -process part can both be estimated, and from them the chronological 



solutions can be displayed. Readers can turn to |Clayton ( 1988) for these solutions 



The cosmoradiogenic chronologies seem to indicate an older galaxy than do the 
direct-remainder chronologies (Clayton] |1988] l. But they also are compromised by 
unique and interesting uncertainties. For the 187 Re cosmoradiogenic chronology, the 
main uncertainty is how greatly its decay rate is speeded by its incorporation into 



stars, where ionization increases its beta decay rate. Yokoi et al. (1983) evaluated 
this effect within an ambitious galactic chemical evolution model that enabled them 
to take into account the fraction of 187 Re that is incorporated into stars and ejected 
again without nuclear processing and how much time such interstellar 187 Re spends 
inside of stars. Their results can be reproduced by a 40% increase of the neutral 
187 Re decay rate. For the Pb cosmoradiogenic chronology, the main uncertainty is 
an especially interesting set of nucleosynthesis problems associated with both the s- 
process part and the r-process part of their abundances. Uncertainty about the Pb/U 
elemental abundance ratio also suggests some caution. 

Certainly it can be said that these chronologies present intricate problems in the 
astronomy with radioactivity. 



Radioactivity and Galactic Chemical Evolution 

In astronomy, the abundances of radioactive species are interpreted through the 
lens of galactic abundance evolution, frequently called galactic chemical evolution 
(GCE). Interpretations of nuclear cosmochronology or of initial solar abundances of 
extinct radioactivities depend on the mean expected abundances in the ISM when 
solar birth occurred. Galactic Chemical Evolution traditionally concerns the chem- 
ical composition of the mean ISM. It features continuous nucleosynthesis between 
a starting time and the time when solar formation occurred. During that lengthy 
galactic period radioactive abundances undergo decay and are incorporated into new 
stars, but they are also replenished by the injection of freshly synthesized radioac- 
tivity and are diluted by low-metallicity matter falling onto the disk. Competition 
among these terms renders simple galactic chemical evolution a showplace for the 
behavior of radioactivity. One differential equation describing that abundance of 
those specific galactic atoms (the red atoms) that will later be incorporated into the 
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solar system can be thought of as the solution of 

dN _ <//>'.;/: N{t) 
dt dt x 



(2.14) 



where T=l/A is the mean radioactive lifetime of the nucleus in question, N(t) is its 
time-dependent abundance within those (red) atoms that are destined for inclusion 
in the solar system, and Z? fl (f ) is the birthrate spectrum of those solar nuclei, defined 
such that dB ^*' dt is the number of stable solar nuclei that were born between times t 



and t+dt. The birthrate spectrum of solar nuclei is unknown, however, so eq. 2. 14 



not be solved. Despite this limitation, workers persisted in attempting to fix the time 



of the beginning of nucleosynthesis from eq 2.9 which students can show does solve 



the differential equation 2. 14 above. Many tried this approach by assuming an easily 



integrable form for dB ^ , most often exp(-Gt), where G is taken (erroneously) to 
be an unknown positive number. That assumption proved seductive to many, because 
the rate of galactic nucleosynthesis is believed to decrease smoothly with time, so 
maybe dB j. ' can be assumed to decline as well. Importantly, however, dB % differs 
greatly from the galactic production rate and actually grows with time. The reason 
for this surprise requires understanding. 

The birthrate spectrum imagines that the nuclei destined for the solar system 
could be tagged at birth (imagine red paint) and ignores all other ISM nuclei. Be- 
cause B & (t) is unknown, however, astrophysicists instead calculate the evolution of 
the interstellar abundances. Their mean values in a well mixed ISM at the time of so- 
lar birth provide the expected solar abundances. The conceptual difference between 
B (f) and the galactic nucleosynthesis rate has confused many unwary researchers 
Many published papers have used their chosen form for dB % and integrated eq 



2.9 



in the hope of fixing the beginning time of galactic nucleosynthesis (the lower limit) 
by comparing those calculated results with known initial radioactive abundances 



in the solar system. A common form because eq 2.9 is then easily integrated is 
= exp(-Gt). In retrospect, however, such poor choices for the form of dB ®^ 
render their conclusions invalid. 

How does GCE clarify the problem of galactic radioactivity? The mass Mq of 
ISM gas is reduced owing to the rate \j/(t) of its incorporation into new stars, but it 
is increased by the rate E(t) at which mass is ejected from old stars and by the rate 
f(t) at which new mass falls onto the galactic disk: 

d ^=- ¥ (t)+E(t)+f(t) (2.15) 

If one takes the ejecta rate E(t) from the spectrum of newly born stars to be a 
fixed return fraction R of the rate \j/(t ) at which mass joins new stars, and if one 
assumes linear models in which the rate of star formation is proportional to the 
mass of gas Mq , the equations governing the interstellar composition of both stable 
and radioactive nuclei can be solved analytically within families of choices for the 



infall rate /(f), as Clayton (1988 1985 i has shown. These analytic solutions for 



24 



2 The Role of Radioactive Isotopes in Astrophysics 



ISM abundances do indeed clarify radioactive abundances. For that purpose Clayton 
writes 

dMr 

-^ = -y{t){\-R) +f(t) = -03M G +f(t) (2.16) 

where a>Mc(t) = (1 — R)y/(t) is valid for linear models. The constant co is the rate 
of consumption of gas by star formation when compensated by gaseous return from 
stars. These are called linear models because the star formation rate \j/(t) is taken 
to be proportional to the mass M G {t) of interstellar gas. That linear assumption is 
not strictly true; but it is plausible in taking the star formation rate to increase as 
the mass of gas increases, and to decline as it declines. Furthermore, it is supported 
by observations of star formation rates in spiral galaxies. From their observations, 
|Gao et aL"] ( |2004| l state, "The global star formation rate is linearly proportional to the 
mass of dense molecular gas in spiral galaxies." The purpose of models of galactic 
chemical evolution insofar as galactic radioactivity is concerned is to understand 
the mean expectation for radioactivity in the ISM. The families of linear analytic 
models are constructed for that purpose. 

Coupled with eq. |2.16| is an equation for the rate of increase of the concentration 
Z of each interstellar nucleosynthesis product. The concentration is defined as the 
mass of species Z in the ISM divided by the total mass Mq of ISM gas and dust. The 
concentration Z, rather than total numbers of atoms, is the quantity traditionally used 
in chemical evolution studies because it is concentration that astronomers measure. 
The mass mz = ZMq of interstellar species Z is governed by 

-^ = -Z\i/(t)+Z E E{t)+Z f f(t) (2.17) 

where Ze and Zf are, respectively, the concentration of Z in the spectrum of stellar 
ejecta, where it is large, and in the infalling gas where it is small. The metallicity in 
inf ailing gas Zf may probably be neglected with good accuracy. 

Many workers have shown that it is a good approximation for analytic under- 
standing to assume that the ejecta E(t) is returned at once from the entire spectrum 
of newly born stars despite the dependence of stellar lifetime on stellar mass. That 
assumption is called the instantaneous recycling approximation, and it is reasonable 
except late in the life of a galaxy when the gas mass and the star-formation rate have 
both become very small. But for galaxies in early and middle lifetime, eq. |2.17| can 
then be written after some straightforward algebra as 

^yc-^-AZ (2.18) 
dt y M G (t) 

where mz — ZMq, Z being the concentration in the ISM gas (taken to be well mixed), 
and y is the yield of element Z, defined as the mass mz of new Z ejected from the en- 
tire spectrum of newly born stars divided by the mass of stellar remnants left behind 
by that entire spectrum of stars. The yield y for primary nucleosynthesis products 
may be taken as a constant despite having small variations in full numerical models. 



The product y(0 in eq 2.18 can be called the galactic nucleosynthesis rate, which, 
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through (0, depends explicitly on the star formation rate {coMc = (1 — R)\jr(t)). 



Equation 2.18 for the ISM concentration differs from eq. 2.14 for the history of 
solar nuclei by the existence of the second term of eq. 2.18 That term reflects the 



loss of interstellar metal concentration Z when ISM gas containing Z collapses into 
new stars and is simultaneously diluted by metal-poor infall. 

Comparison of eq|2.18|with eq|2.14|shows clearly why it is an error to integrate 



eq. 2.14 thinking that the solar-nuclei birthrate spectrum can be mentally equated 



with the spectral rate of galactic nucleosynthesis. Equation 2.18 similarly has the 



galactic nucleosynthesis rate as the first term and the radioactive decay rate as the 



last term; but the astrophysical eq. 2.18 contains the middle term involving both the 



galactic infall rate and the mass of ISM gas, both of which are time dependent. In- 



tegrations of eq 2.14 by assuming B @ (t) ignore the effects of the ISM and can be 
correct only if the assumed form dB Q (t)/dt actually resembles the number of solar 
nuclei born per unit time rather than the galactic nucleosynthesis rate. As the birth- 
date of solar nuclei is usually what is being sought in nuclear cosmochronology, 



one cannot get the answer from eq. 2.14 without first knowing and inserting the an- 
swer. For that reason nuclear cosmochronology must instead be investigated within 
the context of galactic chemical evolution. These points are central to the subject of 
astronomy with radioactivity within solar-system nuclei. 

Clayton ( 1988 ) showed one way forward. One can integrate eq 2.18 analytically 



with the aid of flexible families of functions that enable analytic integration. This 
can be accomplished by form-fitting the ratio f(t)/M(j(t) to an integrable family of 
functions d9 jdt. Such form fitting is much more flexible than it would at first seem 
to be, because most physically plausible time-dependent behavior can be approxi- 
mated by a specific form fitting choice. One such useful choice has been called the 
Clayton Standard Model. It takes 



dd 

dt 



■f(t)/M G (t) 



(t + A) 



(2.19) 



and its explicit functions are given in Appendix Al of this book. 

For the astronomy of radioactivity, the great merit of this approach is exact ana- 
lytic solutions for Z(f ) for both stable and radioactive isotopes while simultaneously 
yielding exact functional representations of the mass M(j(t) of interstellar gas and 
of the infall rate /(?). Understanding these analytic models greatly aids understand- 
ing the behavior of radioactivity within more general numerical models of galactic 
chemical evolution. The numerical approach is to place the evolution of the galaxy 
on a computer, taking into account the evolutionary lifetime of each star formed and 



the specific nucleosynthesis products to be ejected from each star (e.g. Timmes et al. 



1995 ). Although this approach is undoubtedly correct, it obscures theoretical under- 
standing that can be seen more easily within analytic models. Furthermore, surveys 
of nuclear cosmochronology can more easily be carried out within analytic models. 

All well mixed models of galactic chemical evolution can at best yield only an 
average expectation for the ISM. The true ISM is inhomogeneous in space and the 
nucleosynthesis rate is sporadic in time rather than maintaining its steady average. 
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Nonetheless, well mixed models, both analytic and numerical, are important in lay- 
ing out the results that would be true for a rapidly mixed ISM and a smooth rate of 
galactic nucleosynthesis. One interprets the observations of radioactive abundances 
against the backdrop of that expectation. 

Another result of great importance for short-lived galactic radioactivity is best 
illustrated within Clayton's standard model. For short-lived radioactivity, i.e. when- 
ever X is small, the mean concentration in the galactic ISM is 



Zl =yco[X + — j (2.20) 

From this equation and one for a stable isotope the abundance ratio of a short-lived 
radioactive nucleus to that of a stable primary isotope in the ISM is larger by the 
factor (k + 1) than one might estimate without taking galactic chemical evolution 
into account. As an example, consider the interstellar ratio of radioactive 26 Al to 
stable 27 Al. The formula 

Z^-yW( k+1 )?* (221) 
Z( 27 A1)- t e (Z21) 

1993). This result is larger by the 



can be derived from such a ratio (Clayton et al. 



2.10 



which had traditionally been used 



factor (k+ 1 ) than an estimate using only Eq. 
in oversimplified discussions of the steady-state amount of short-lived interstellar 
radioactivity. That oversimplification can be found in almost all published papers 
on short-lived interstellar radioactivity. The extra factor (k+ 1), which is an effect 
of infall of low-metallicity gas, became important when a large mass of interstellar 
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Al was detected by gamma-ray astronomy (see Chapter 7). 



2.2.2 Gamma-Ray Lines from Galactic Radioactivity 

One thinks of at least three reasons that the idea of astronomically detecting galac- 
tic radioactivity did not occur until the 1960s. In the first place, MeV range gamma 
rays do not penetrate to the ground, but must be detected above the atmosphere. 
Secondly, detecting MeV gamma rays and measuring their energies is a quite dif- 
ficult technology, especially so when the background of cosmic-ray induced events 
above the atmosphere is so large. Thirdly, even a back-of-the-envelope estimate of 
expected rate of arrival from a stellar source is discouragingly small. 



The Rice University Program 

In 1964 a new aspect of astronomy with radioactivity arose. Robert C. Haymes 
was hired by Rice University, and he spoke there with Donald Clayton about the 
possibility of sources of galactic radioactivity that he might seek with an active an- 



2.2 Disciplines of Astronomy with Radioactivity 



27 



ticoincidence collimation for a Nal detector flown beneath a high-altitude balloon. 



Burbidge et al. ( 1957 1 had speculated (incorrectly) that the exponential 55-night de- 



cline of the luminosity of many Type la supernovae was the optical manifestation 
of the decline of the spontaneous-fission radioactive decay of 254 Cf in the ejecta of 
the supernova, which would quickly have become cold without some heating mech- 
anism that had to decline with a roughly two-month halflife. The 254 Cf nucleus 
would be synthesized by the r process (assumed to occur in Type la supernovae) 
along with uranium and thorium. The large kinetic energy of its spontaneous-fission 
fragments would be converted to optical emission by being degraded by atomic 
collisions within the ejected gas. This was called the californium hypothesis. But 
B 2 FH had said nothing about gamma-ray lines. The novel excitement at Rice Uni- 
versity was the new idea that concepts of nucleosynthesis could be tested directly if 
the associated gamma-ray lines from the r-process radioactivity could be detected 
on earth coming from supernova remnants. Haymes estimated rather optimistically 
that his detector could resolve lines having flux at earth greater than about 10~ 4 



cm 2 s 1 . The first scientific paper written with that goal was soon published ([Clay- 



ton and Craddock||1965| l. It evaluated the full spectrum of radioactivity by an r pro- 
cess normalized to the yield proposed by the Cf hypothesis. They found several 
promising lines. For example, the strongest from a 900-year old Crab Nebula would 
be 249 Cf, presenting 10 _4 cm _2 s _1 gamma-ray lines having energy 0.39 MeV. It is 
no accident that the 249 Cf halflife (35 1 yr) is of the same order as the age of the Crab 
Nebula. Given a broad range of halflives in an ensemble of nuclei, the one giving 
the largest rate of decay has mean life comparable to the age of the ensemble. It 
can be confirmed that in a remnant of age T the rate of decay per initial nucleus 
is XzT^ 1 ', which is easily shown to be maximal for that nucleus whose decay rate 
A = 1/t = l/7\ A photograph (Fig. 2.1) from those Rice years is included. 

Realistic expectations were not so sanguine, however. The r-process yield re- 
quired by the Cf hypothesis appeared excessively large (Cla yton and Craddock| 
1965[ >. The r-process nuclei in solar abundances would be overproduced by a factor 



100 if all Type la supernovae produced that yield. Despite that reservation, an ex- 
citing chord had been struck. A new astronomy of radioactivity appeared possible, 
one having significant scientific payoff, and Haymes's gamma-ray telescope began 
a program of balloon flights. The program did not discover r-process radioactiv- 
ity, but it did discover apparent positron-annihilation radiation from the center of 
our galaxy ( |Johnson et al. 1972| l and hard x-rays from several sources. The obser- 
vational program at Rice University also trained two of the principal investigators 
(G.J. Fishman and J.D. Kurfess) and two of the project scientists (W.N. Johnson and 
C. Meegan) on NASA's later Compton Gamma Ray Observatory, which went into 
orbit in 1991. And Clayton began a NASA-sponsored research program for laying 
out theoretical expectations for nucleosynthesis produced radioactivity and its as- 
sociated nuclear gamma-ray lines in astronomy. Gamma-ray-line astronomy was in 
that sense launched by the Rice University program. 

At least two other groups followed the Rice University lead and became and 
remained leaders for the subsequent decades. They constructed differing telescope 
techniques from that of Haymes Nal scintillator. One was the JPL research group 
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Fig. 2.1 Robert C. Haymes and Donald D. Clayton with Haymes' gondola for his balloon-borne 
gamma ray telescope at Rice University in 1973. Haymes pioneered gamma-ray astronomy in the 
MeV region and Clayton developed gamma-ray-line targets for this goal and for nucleosynthesis. 
This program at Rice University trained two Principal Investigators and two Project Scientists for 
later experiments on NASA's Compton Gamma Ray Observatory. 



- A.S. Jacobson, J. C. Ling, W. A. Wheaton, and W.A. Mahoney. Using a cluster 
of four cooled high-purity germanium detectors on NASA's HEAO 3, which was 
launched in 1979, they discovered the first galactic gamma-ray line from radioac- 
tive decay. A group at MPI fur Extraterrestrische Physik in Munich, V. Schonfelder, 
P. von Ballmoos and R. Diehl, developed a Compton-scattering telescope capable 
of good energy resolution and a greater angular view of the sky. Their active results 
from balloon-borne launches in the 1980s presaged the splendid results of their in- 
strument on NASA's Compton Gamma Ray Observatory. Jacobson's JPL group was 
not so lucky, as their experiment was removed from CGRO as a cost saving. Other 
groups joined the effort to create this new astronomy with radioactivity. 

Earliest Predictions of Detectable Gamma-Ray Lines 

Inspired by the Rice University balloon program, a series of studies of the nucle- 
osynthesis of radioactivity was undertaken. It may seem surprising in retrospect 
that this could not have been done by a routine search of the chart of nuclides in 
conjunction with nucleosynthesis theory. But the data base necessary did not ex- 
ist. Four decades ago, nucleosynthesis theory was undergoing rapid expansion and 
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clarification, so that it was a series of new insights that laid out each new prospect 
as it was discovered. Unquestionably the most important of these developments oc- 
curred in 1967 when a Caltech group (Bodansky et al. 1968a|b i and a Yale group 
(Truran et al. 1967 1 showed with differing research techniques that the e process 
of Hoyle et al. ( 1 946 1 and Burbidge et al. ( 1957 1 for the nucleosynthesis of the 
iron abundance peak was incorrectly applied. Instead of producing iron isotopes 
as themselves within a neutron-rich nuclear-statistical equilibrium (NSE), as Hoyle 
had suggested, they were produced as isobars of radioactive nickel that decayed to 
iron isobars only after nucleosynthesis and ejection from supernovae. Owing to the 
increasing strength of the Coulomb energy in nuclei, the most stable nucleus for An 
nuclei shifts from Z = N below A=40 to N = Z + 4 for A=44-60. But insufficient 
neutron excess required those major nuclei to form as unstable Z = N nuclei. This 
brought into view the astronomical significance of several important radioactive nu- 
clei. A historical study, Radiogenic Iron, was written later to lay out the many sig- 



nificant astrophysical issues that hinged on this new understanding ( Clayton et al. 
[1993) . 

It is of interest in a historical sense to consider the sequence of subsequent 
predictions of target radioactivity for astronomy. Clayton et al. ( 1969) jumped on 
the newly discovered radioactivity in the revised equilibrium processes (preceeding 
paragraph). Easily the most important nucleus was 56 Ni, having Z = A, whose beta 
decays to 56 Co and thereafter to 56 Fe, having A = Z+4, are accompanied by numer- 
ous gamma ray lines that might be resolved by spectroscopic gamma-ray detectors. 
G. J. Fishman was a research student in R. C. Haymes' Rice University balloon pro- 
gram when he and Clayton began calculating this spectrum of lines and their time 
dependences. Their motivation was to test nucleosynthesis theory, especially the 
fresh new view of explosive nucleosynthesis that had led to 56 Ni nucleosynthesis. 
S. A. Colgate joined this effort because he had independently begun to investigate 
whether the energy deposition from the radioactivity could explain the exponen- 



tially declining Type la supernova light curves (Colg ate et al.| [1969). Burbidge et al. 
( |1957| l had attributed those long lasting light curves to heating by s pontaneous fis- 
sion fragments from freshly produced 254 Cf. As the simplest model ( Clayton et al 



1969| > evaluated fluxes at earth from the rapid expansion of a solar mass of concrete, 
much of which had been explosively processed to 56 Ni in the way anticipated by 



Bodansky et al. ( 1968a b I. The lifetime of 56 Co was seen as favorable for sufficient 



expansion for the gamma rays to escape. They neglected a massive envelope, as 
might overlie Type II cores, realiz ing that such an env elope would seriously limit 
the escape of the 56 Ni gamma rays. Clayton et al. ( 1969 1 was a very important paper, 
not only because it was correct in its prediction of the strongest gamma-ray nuclear 
source, but also because of its galvanizing effect upon experimental teams and on 
NASA planning. It later was included in the AAS Centennial Volume of seminal 
papers of the 20 th century. NASA Headquarters suggested funding of a research 
program at Rice to lay out additional prospects for this new astronomy of radioac- 
tivity. The NASA grant at Rice was entitled Prospects for Nuclear Gamma-Ray 
Astronomy. But the detection of those 56 Co gamma rays did not become possible 
until SN1987A occurred in 1987. 
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|Clayton et al.| ([T969) also first suggested that unknown galactic supernovae (un- 
seen owing to optical obscuration) may be discovered in gamma-ray lines. The lines 
from radioactive 44 Ti, with 60-yr halflife, ought to be detectable from several super- 
novae that have exploded during the past two centuries, considering that about three 
per century seem to occur on average within the Milky Way. This exciting idea has 
great implications for galactic nucleosynthesis and astrophysics; but its hope has 
been frustrated for astrophysical reasons that are not yet understood. 

It was quickly realized (Clayton and Silk 1969) that the entire universe might 
be filled with detectable 56 Co gamma rays from the past history of supernova ex- 
plosions throughout the universe. This work suggested that observable gamma rays 
owing to the collective effects of many supernovae rather than specific young super- 
novae might be targeted. Detecting the universe in this way remains a goal of the 
astronomy of radioactivity. 

Two years later it was proposed (Clayton 1971 j > that the collective effects of 
galactic supernovae might allow one to observe gamma-ray lines from long-lived 
radioactive nuclei whose emission would be too weak from individual supernovae. 
That first work in that regard focused on 60 Fe, whose 2.6 Myr halflif^jmakes emis- 
sion from an individual supernova too slow to be observable for the modest number 
of 60 Fe nuclei produced within a single supernova. The 60 Fe nucleus emits a 59 ke V 
gamma ray upon decay, and its daughter 60 Co emits gamma-ray lines of 1.17 and 
1.33 MeV. Reasoning that during its long mean lifetime some 50,000 supernovae 
occur in the Milky Way, their collective effect should be observable. This reasoning 
applied equally well thirteen years later to the first interstellar radioactivity to be 
detected, that of 26 Al. 

Surprisingly in retrospect, several years passed before [Clayton ( 1974) realized 
that when 57 Co, the daughter of 36-hr 57 Ni, decays, it also emits favorable gamma- 
ray lines and that its longer 272-day halflife would cause it to radiate these after the 
56 Co was almost extinct. This made 57 Co much more significant than had been ap- 
preciated. Clayton ( 1974 1 drew attention to the significance of 57 Co and to several 
new ideas for astronomy with radioactivity. Firstly, given an appropriate galactic su- 
pernova one might measure the material thickness overlying the radioactive cobalt 
by the distinct times at which differing gamma-ray lines reach maximum flux. This 
happens because the structural opacity for the 2.60 MeV gamma ray is only about 
half that for the 0.84 MeV gamma ray (both from 56 Co decay), so that the 2.60 MeV 
gamma ray should peak earlier despite the equal rates of emission. Additionally, the 
0.12 MeV gamma ray from the slower 57 Co decay suffers even more absorption and 
so peaks later yet in time. These time delays have not yet been measurable because 
a time-dependent gamma line flux requires a bright supernova which would have 
to occur within our own galaxy. Clayton ( 1974| l also warned that hydrodynamic in- 
stabilities could allow the central material to burst through in streams and thereby 
appear earlier than anticipated. Early escape of lines from 56 Co was later detected 
in SN1987A, which exploded in the Magellanic Clouds, and from which gamma- 
ray lines from both 56 Co and 57 Co were first detected (see Chapter 3). Earliest of 



2 At that time, the halflife of 60 Fe was best-known as 1.5 My 
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the recorded 56 Co photons were by the gamma ray spectrometer aboard NASA's 
Solar Maximum Mission, which happened to be in orbit when SN1987A exploded 
and whose sun-pointing spectrometer was reinterpreted as a gamma ray telescope 
(Leising et al. 1990|). Several balloon-borne gamma-ray spectromete rs were flown 



and also detected sb Co lines. When 57 Co was detected in SN1987A by Kurfess et al. 



( |1992| l with the Oriented Scintillation Spectrometer Experiment on NASA's Comp- 
ton Gamma Ray Observatory during summer 1991, four years (about six halflives 
of 57 Co) had elapsed since the explosion, so that the 57 Co abundance had decayed 
to only about 1.5% of its initial amount. Tension surrounded the hope that it would 
still be detectable. It was, with 4a significance, implying that the initial abundance 
ratio 57 Ni/ 56 Ni was about twice the ratio 57 Fe/ 56 Fe measured in terrestria l iron . 

1974} pro- 



At almost the same time as the 57 Co prediction, |Clayton and Hoyle 



posed that gamma rays may be detectable from the more common nova explosions 
owing to the radioactivity created by the nova outburst. These involve the positron- 



annihilation line from the hot-CNO burning that powers the nova outburst (Star- 



rfield et al.] |1972| l. These 511 keV photons would have to be detected very quickly, 
within roughly 10 3 s after the outburst because of the short halflives of CNO radioac- 
tive nuclei. They also pointed out the prospect of detecting a 1 .274 MeV gamma ray 
following 22 Na decay. Because of its longer halflife (2.6 yr), 22 Na remains alive for 
the roughly 10 6 s required for the nova ejecta to become transparent to the gamma 
ray. These goals have not succeeded yet, but they remain a realistic hope of measur- 
ing the thermonuclear power of the nova. The model of the nova makes it a remark- 
able laboratory for thermonuclear explosions. Detection could confirm the model or 
rule it out. Some years later the possibility arose that the 478 keV gamma-ray line 
following the decay of radioactive 7 Be to 7 Li might also be detectable from novae 
(Clayton 1981 ). Detectabilty requires the nova envelope accreted from the compan- 
ion star to be enriched by an order of magnitude in 3 He, but that is perhaps to be 
expected. If so, the nuclear reaction 3 He + 4 He =>• 7 Be + y creates the radioactive 
nucleus in the thermonuclear flash. One positive aspect is that the 7 Be halflife is 
large enough for the envelope to become transparent. 

The radioactive 22 Na nucleus is also produced in supernovae, although the com- 
plicated details of its nucleosynthesis therein were not well understood when |Clay~^| 
|ton| ( |1975c] l proposed its detectability within supernovae. It might be observable for 
a decade, allowing ample time to become transparent to the 22 Na gamma-ray line. 

Almost the last good prospect to be predicted turned out to be the first actually 
observed! That was the radioactive 26 Al nucleus (Ramaty and Lingenfelter 19771. 
Like radioactive 60 Fe, the 26 Al nucleus is long-lived and thus decays too slowly to 
be detectable from individual explosions. It is the cumulative yield of many events 
over the past Myr or so that was discovered. This is the topic of the next section. 

The predictions whose history has been recounted here placed tantalizing targets 
before the community of experimental physicists. Many of these assembled teams 
to design, build and fly gamma-ray spectrometers capable of detecting these lines. 
NASA entertained a grants program pursuing these goals, which create a new wave- 
length range for astronomy. Moreover, the will to fund the Compton Gamma Ray 
Observatory was strengthened by these specific hopes. That hope was fueled also 
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by the surprising discovery of a gamma-ray line from 26 Al nuclei in the interstellar 
medium. 



The Surprise 26 Al 



The first detection of an interstellar radioactive nucleus came as a surprise. The his- 
tory of this radioactive isotope in astrophysics reveals that correct predictions can 
be made on the basis of inadequate reasoning, that scientists may be blindsided by 
their own excessive faith in their own pictures and beliefs, and that experimental 
discovery is the arbiter. The possibility of detecting 26 Al nuclei by observing the 
1.809 MeV gamma-ray line emitted following its decay to 26 Mg was suggested by 
Ramaty and Lingenfelter] ( | 1 977| > and independently by |Arnett et al.| ( |l977| ). Their 



interesting papers were not quantitative predictions because the grounds for their 
suggestions were not convincing. They suggested wrongly that interstellar 26 Al nu- 
clei would be detectable if the production ratio in massive stars is P( 26 Al)/P( 26 Mg) 
= 10~ 3 . However, that production ratio, which was expected by carbon burning in 
massive stars, was inadequate for detectability. The argument showing that that pro- 
duction ratio was unobservable ( Clayt on] |1984]l was submitted for p ublication only 

4J|. The predicted 



after interstellar 26 Al was discovered (Mahoney et al. 



1982 



1984) 



1.809 MeV flux was too small by a wide margin to be detected by contemporary 
gamma-ray spectrometers. So something was wrong. 

Using the gamma-ray spectrometer on NASA's third High Energy Astrophysi- 



cal Observatory (HEAO 3) |Mahoney et al. ( 1982 1984 1 reported a measured flux 
of 1.809 gamma rays that required about 3M of interstellar 26 Al nuclei. That 



quantity was far larger than the predicted interstellar mass. Nonetheless, both Ra- 
maty and Lingenfelter| ( |1977 1 and Arnett et al. ( 1977 1 had urged, on general grounds 



and on intuitive arguments, that the 1.809 MeV interstellar gamma-ray line be 
sought. It was indeed observable, the first to be detected. 

What were the conflicts that this discovery illuminated? Using Eq. 
number of radioactive nuclei born during the last mean lifetime (T= 
of 26 Al nuclei, the time-average interstellar abundance ratio for aluminum isotopes 



ZTOjforthe 
1.05xl0 6 yr) 



was traditionally estimated to be 2b Al/ 2/ Al = f( 26 Al)/P( 2/ Al) r/10 lu yr = 10"' if 
one uses f( 26 Al)/P( 27 Al)=10 -3 . Since total ISM mass is about 10 10 M o , it contains 
about 5.8 x 10 5 M o of stable 27 Al. Multiplying by the above isotope ratio, the ISM 
would carry about 0.06 M of 26 Al nuclei, woefull y inadequate for the observed 
3M of interstel lar 2 6 Al (|Mahoney et al.[[l982]|l984j >. So although Ramaty and Lin- 
genfelter ( 1977 1 and |Arnett et al. ( 1977 1 correctly suggested that Al nuclei might 



provide a suitable ISM radioactivity to seek, their reasoning did not convincingly 
justify that hope. The discovery of so much interstellar 26 Al was surprising and 
meant that some assumptions were quite in error. In fact, several were in error. 

Confirming detections of the 26 Al radioactivity were quickly made from balloon- 
borne gamma-ray spectrometers ( |von Ballmoos et~aL 1987| MacCallum et al. 
1987 ), and the total mass was later measured accurately by the gamma ray spec- 
trometer aboard NASA's Solar Maximum Mission, whose sun-pointing spectrometer 
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had been reinterpreted as a gamma ray telescope (Share et al. , 1985 ). The existence 
of 2-3 M , of 26 Al was then beyond doubt. The 1991 launch of NASA's Comp- 
ton Gamma Ray Observatory with its imaging COMPTEL Compton telescope (see 
Chap. 10) produced the m ost detailed data about the spatial distribution of interstel- 
lar 26 Al |Diemetd][l995 1. 

Stepping back to the history of the mid 1970s, a related issue for the astronomy 
of 26 Al radioactivity had arisen about one year prior to Ramaty & Lingenfelter's 
1977 suggestion; namely, it was discovered (Gray et al. 1974 Lee et al. 1977] > 



that the molecular cloud from which the sun had formed apparently carried within 
it the large isotopic ratio 26 Al/ 27 Al=5x 10~ 5 at the time the planetary system be- 
gan to form. Had that ratio been applicable to the ISM as a whole, it would have 
corresponded to 29 M of radioactive 26 Al nuclei-very much more than was dis- 
covered in 1982 by Mahoney et al. The seeming impossibility of producing such a 



large quantity by nucleosynthesis prompted the idea (Cameron and Truran 1977 



Wasserburg et al. 1982 ) that a supernova within the molecular cloud from which 



the sun was born had injected the solar 26 Al radioactive nuclei into the collapsing 
solar portion of the cloud just prior to the solar birth. In that case the solar cloud was 
quite atypical of the ISM at large. The initial solar 26 Al nuclei in meteorites became 
the standard bearer for a class of extinct radioactive nuclei that required special lo- 
cal production realted to solar birth. But owing to the huge consequent ratio in the 
forming sun, it seemed to bear no transparent relationship to the existence of several 
solar masses of interstellar 26 Al nuclei, which necessarily represents nucleosynthe- 
sis of 26 Al nuclei by many thousands of supemovae spread out in time over 1-2 Myr 
prior to today. This discovery from meteorites brought the idea of interstellar 26 Al 
nuclei to the attention of astrophysicists. 

What new ideas brought the expected 0.06M Q of 26 Al nuclei in the ISM into 
line with the observed 3M ra of 26 Al nuclei? Some suggested that novae ( Clayton 



1984, Woosley 'and Weaver] 1 1980) or AGB stars, rather than supernovae, were the 
source of the 2b Al nuclei; but that possibility was ruled out by observations made 
later with the Compton Telescope (COMPTEL) following the launch of Compton 
Gamma Ray Observatory. Those observations (Die hl et al.| 1995[ ) showed spatial 
concentrations of 26 Al nuclei in ISM regions where star formation is currently ac- 
tive. Spatial correlation of 26 Al nuclei with massive stars was clear. Attention there- 
fore returned to supernovae sources with the realization ( |Arnould et al.| |1980[ l that 
hot hydrogen burning shells of massive stars constitutes a significant extra source of 
26 Al. The small production ratio in carbon burning (P( 26 A1)/P( 27 A1) =10~ 3 ) must 
be augmented by 26 Al-rich convective shells in massive stars. Those calculations 
became a major industry with numerical models of presupernova evolution. Be- 



cause Al/ Al ratios of order 0.1 exist in these shells, their contribution calcu- 
lated by using time-dependent numerical models of presupernova evolution would 
be needed. Weaver and Woosley ( 1993 i calculated that the average production ra- 
tio from a standard spectrum of massive stars was near P( 26 A1)/P( 27 A1)=0.006, six 
times larger than the earlier estimates. Secondly, the estimated ratio of mean inter- 
stellar 26 A1/ 27 A1 must be increased because most of the stable 27 Al is locked up in- 
side old stars whereas the live 26 Al nuclei will still be overwhelmingly in the gaseous 
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ISM. This reasoning augments the expected ratio by the factor (k+1) derived from 



the standard analytic models of chemical evolution of the mean galaxy. Clayton et al. 
( 1993 ) presented that argument for the expected interstellar ( 26 A1/ 27 A1) abundance 
ratio. Because k=2-4 seems likely from other astrophysical arguments, this (k+1) 
correction also increases the interstellar ratio by another factor 3-5. These two ef- 
fects combined amounted to an increase of a factor 20-30 in the initial expectations, 
bringing the expected value from 0.06 M & to 1-2 M Q of 26 Al nuclei in the ISM. The 
original conflict with global theory was largely resolved. 

Moreover, the largest observed flux from 26 Al nuclei concentrated in ISM re- 
gions where star formation is currently active. This directly challenged the theoreti- 
cal simplification that adopts a star-formation rate that is constant in time and spread 
smoothly through the spatial ISM. Models of chemical evolution of the galaxy use 
that simplification in order to be calculable. But the observations of 1.809 MeV 
gamma rays painted a picture of star formation occurring sporadically whenever 
and wherever large regions of ISM become vulnerable to prolific star formation. 
The massive 26 Al emission regions stunned and surprised many experts (including 
this writer). The fluxes from these are moderated by the inverse square of their re- 
spective distances, further complicating discussions of the mass of interstellar 26 Al. 
The history of 26 Al nuclei in the ISM provides a textbook example of scientific 
progress into a new discipline. Experts disagree in their predictions; correct predic- 
tions often rely on intuition as much as on scientific justification; surprises often 
attend exploratory experimental surveys; and scientists from many disciplines am- 
plify the relevant ideas and measurements. Astronomy with radioactivity is such a 
science, and its first discovery of an interstellar gamma-ray line emitter, the surprise 
of 26 Al, was bathed in confusion. 



2.2.3 Radiogenic Luminosity 

A big problem lay in the path of trying to understand how explosive objects in as- 
tronomy could remain bright. Their luminosity was expected to fall rapidly as the 
objects expand. Suppose a rapidly expanding object consists of an ideal nondegen- 
erate gas and that the expansion is too rapid for gain or loss of heat; that is, the 
expansion is adiabatic. The internal heat of such an object is rapidly lost to the me- 
chanical work of the expansion. For adiabatic expansion TV^ 1 = constant during 
the expansion, where y is the ratio of specific heats at constant pressure cp to the 
specific heat at constant volume cy- For a perfect monatomic gas y= cpjcy = 5/3. 
Thus TV 2 / 3 = constant for such an expanding hot object. If the expanding object 
can escape gravitational binding, the radius of an exploding object may be approx- 
imately R = vt, so that volume V is proportional to ? 3 . Thus one would expect that 
Tt 2 is constant during the expansion. The problem is then that if luminosity L is 
proportional to R 2 T 4 , as in a black body, one expects L to be proportional to f~ 6 . 
However, many supernovae that are bright after one week are still bright after three 
weeks instead of dimmer by the large expected factor near 3~ 6 . Similar problems ex- 



2.2 Disciplines of Astronomy with Radioactivity 



35 



ist for novae expansions. The problem, then, lay in discerning what source of energy 
could keep the expanding objects hot. The answer is radioactivity. The radiogenic 
luminosity of exploding objects became one of the major aspects of astronomy with 
radioactivity. 



Radiogenic Luminosity in Supernovae 

It was excellent luck that SN 1987A exploded nearby in the Magellanic Cloud in 
1987. It became the most observed explosive astronomical event of all time. As- 
tronomers recorded its emissions in every possible wavelength band, which proved 
necessary in exposing how radiogenic luminosity worked in that kind of supernova. 
The observations showed that after July 1987 the total power output of SN 1987A 
declined exponentially for about two years, accurately tracking the 77.2-day halflife 
of 56 Co. This confirmed that the energy of the positrons emitted and of the subse- 
quent gamma-ray emission following each 56 Co decay were efficiently converted to 
optical and infrared luminosity. This old idea for declining supernova light c urves 
was first treated for the identification of 56 Co as the relevant radioactivity by Col- 



gate et al. ( 1969| l. The observed value of the total luminosity of SN 1987A required 
O.O75M of 5f> Ni to have been synthesized within the core matter that escaped the 
central neutron star. With each scattering of a gamma ray it loses roughly half of its 
energy to the recoil electron, which is quickly degraded into heat. Once its energy 
has been scattered below 5 keV the large photoelectric opacity converts its entire 
remaining energy to electron energy, which is also degraded into heat. 

Supernovae come in two main types, with additional structural subtypes. Type 
la and Type II differ in how compact they are, how much radioactivity is produced, 
how that radioactivity is distributed, and how much envelope overlies the radioac- 
tivity. The nucleosynthesis of 56 Ni is easily the largest and most important source of 
heating power for the radiating gas. Such issues impact the way in which radioac- 
tivity provides luminosity for the remnant. To truly judge the degree of scientific 
understanding of radiogenic luminosity requires astronomers to infer, to the best of 
their ability, the structure of the exploding object. The luminosity being observed 
has complicated relationship to the location of the radioactivity. In SN1987A itself, 
for example, the presupernova star was blue, implying smaller size than the com- 
mon red presupernova stars. The amount of overlying matter was modest owing to a 
large extent of mass loss during the prior evolution of the star. These characteristics 
influenced the relative importance of the radiogenic luminosity to the luminosity 
caused by the shock heated envelope. 

Although 56 Ni provides the largest and most important source of heating power, 
other radioactive nuclei also play very important roles. The 56 Co daughter of 56 Ni 
is actually more observable in its effect on radioactive luminosity. This is primar- 
ily because its 77.2-day halflife allows it to remain alive after the 6.08-day parent 
56 Ni has decayed. Therefore, 56 Co delivers heating power at a later time when the 
expanding remnant is much larger and has suffered more adiabatic cooling. In these 
circumstances, quantitative astronomy with radioactivity requires careful evaluation 
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of how, when and where the energy released by radioactivity is deposited in the 
expanding gas ( The et al. 1990 1. 

The nucleosynthesis of 5 / Ni also plays a role in radiogenic luminosity because 
it decays rather quickly to 272-day 57 Co. Owing to its longer halfiife, 57 Co is still 
providing radioactive power when the 56 Co has decayed to negligible abundance. 
This transition of dominant radioactive power occurs after about two years. When 
Compton Gamma Ray Observatory (Kurfess et al. 1992) measured with OSSE the 
actual mass of 57 Co in SN1987A, it corresponded to twice the solar ratio 57 Fe/ 56 Fe. 
This surprised many who had inferred its abundance to be five times the solar ratio, 
a large isotopic ratio that had been deduced on the basis of the bolometric lumi- 
nosity exceeding t he instantaneous power from 56 Co decay ( Suntzeff et al.[[l992[ 



Dwek et al. 



1992 1. The radiogenic luminosity truly deriving from 57 Co decay was 



demonstrated by the OSSE measurement to have been overestimated by a factor 
near 2.5. That discovery team (Clayton et al. 1992) > then advanced a new theoretical 
aspect of radiogenic luminosity, what they called delayed power. They presented a 
model showing that when the remnant gets sufficiently dilute owing to its expan- 
sion, the rate of electronic recombination cannot keep up with the rate of ionization 
that had been produced somewhat earlier by residual amounts of 56 Co radioactivity. 
Because of the high degree of ionization, the radiogenic luminosity from 56 Co ra- 
dioactivity begins to exceed the instantaneous rate of energy deposition from 56 Co 
decay. In effect, the radiogenic luminosity reflects the rate of 56 Co decay from a 
somewhat earlier time than the time of observation. This observed upturn in lumi- 
nosity had been attributed to 57 Co instead of delayed power from 56 Co decay. Such 
time-dependence became a new aspect of astronomy with radioactivity. 

A more general aspect of radiogenic emission lies in any astronomical investi- 
gation in which the gas possesses a higher degree of energy excitation than could 
be expected in the absence of radioactivity. For example, the hydrogenic lines of 
He + in SN 1987A were interpreted as a consequence of the high-energy photons 
escaping from the supernova interior. Hard X rays emanting from the expanding 
supernova 1987A were another clear example of Compton-scattered radioactivity 
gamma rays escaping from the interior; and only radioactivity can produce so many 
hard X rays. Another example might be the detectable presence of doubly ionized 
species, an ionization state that would make no sense in a purely thermal setting. 
All such radiogenic possibilities are aspects of radiogenic luminosity, and as such a 
part of astronomy with radioactivity. 

A similar transition between identities of powering radioactive nuclei occurs 
when the decay power of 60-yr 44 Ti exceeds that of 57 Co. And no doubt there exists 
a similar problem, namely, does the rate of radiant emission of 57 Co power keep 
up with the rate of 57 Co decay power. The radioactive 44 Ti is created in almost the 
same location as are the Ni isotopes. This occurs in the explosive burning of oxygen 
and silicon (Bodan sky et al.|[l968a|b[|Woosley et al. 1973} and reflects the extent 
to which 28 Si has been decomposed into the silicon-burning quasiequilibrium. But 
radiogenic luminosity from this 44 Ti nucleus depends even more strongly on how 
its emissions are converted to luminosity and the extent of the time lag between ra- 
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dioactive decay and the creation of observed photons. Its contributions to delayed 
radiogenic luminosity require careful assessment of those issues. 

Such considerations are significant in the scientific understanding of astronomy 
with radioactivity. The luminosity of a young supernova remnant is one of the eas- 
iest observations of radioactivity in astronomy; but its interpretation requires fairly 
accurate description of the structure of the exploding object. 



Radiogenic Luminosity in Novae 



The brightening of a nova explosion and the associated ejection of matter from its 
expanded envelope are also issues that are dependent upon the deposition of ra- 
dioactive heat to the envelope (Starrfield et al. 1972). Those authors showed that 
a successful model needs to mix carbon from the surface of the underlying white 
dwarf into the hydrogen-rich accreted envelope in order to have sufficient radioac- 
tivity produced by the thermonuclear flash that triggers the nova event. That burn- 
ing is the hot CN cycle. They showed that radioactive decay keeps the envelope hot 
while it expands to larger radiating surface area and therefore increased luminosity. 
They also showed that without the radioactive afterburner, mass ejection would not 



occur. Truran ( 1982) gives more nuclear-physics details of the nova explosions. 



2.2.4 Extinct Radioactivity and Immediate Presolar 
Nucleosynthesis 

Xenology Revisited 

John Reynolds had discovered in 1959 that extinct radioactive 129 I appeared to have 
existed initially in meteorites. The observed ratio to stable I is near 
This stood for decades as a nearly unique example. Then in the mid-1970s new 
discoveries and ideas greatly enlarged the context of astronomy with radioactiv- 



ity. Using Eq. 2.21 for the number of radioactive nuclei born during the last mean 
lifetime (t= 23.5 x 10 6 yr) of 129 I nuclei, the mean interstellar abundance ratio for 
iodine isotopes is expected to be 

129 I/ 127 I = P( 129 I)/P( 127 I)(k + l)T 129 /t0 = 0.013 -0.022 

using P( 129 I)/P( 127 I) = 1.5, ?0=8 Gy for the presolar duration of nucleosynthe- 
sis, and k = 2 — 4 within the standard model of radioactivity in galactic chemical 
evolution. The measured ratio is very much smaller than this expected ratio. The 
interpretation was that the solar molecular cloud sat dormant, with no new nucle- 
osynthesis input, for 4-6 mean lifetimes of 129 I to allow decay to reduce its activity 
to the measured level. As plausible as this waiting period seemed to be for the next 
two decades, it was not correct. 
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Reynolds' laboratory also discovered that that extinct radioactive 244 Pu appeared 
to have existed initially in meteorites. It was measured by an anomalous pattern of 
several Xe isotopes that was consistent with the spontaneous-fission spectrum from 
244 Pu. (Hudson et al. 1900[ > presented more modern data setting its abundance rel- 
ative to that of 2i2 Th as 244 Pu/ 232 Th = 3 x 10~ 3 . The many studies of these first 
two extinct radioactivities came to be called xenology because of their reliance on 
excess isotopic abundances of Xe isotopes in meteorites. With a production ratio 
P( 244 Pu/ 232 Th) = 0.5 and mean lifetime T 244 = 1 15 X 10 6 yr, the expected inter- 
stellar ratio would be 0.02-0.03 - larger by tenfold than the solar amount. But in 
this case a sequestering of the solar cloud to allow 129 I nuclei to decay would not 
sufficiently reduce the 244 Pu abundance. Trying to resolve the xenology puzzles oc- 
cupied many meteoriticists for two decades. Then three new developments altered 
xenology studies. 

Firstly, existence of such 129 Xe-rich and fission-Xe-rich interstellar dust was pre- 
dicted (Clayton 1975a I by the first argument that dust containing evidence of ex- 
tinct radioactivity should condense within supernovae and influence the early solar 
system. This was an exciting new idea in the astronomy of radioactivity. However, 
it appears that chemical ways of forming meteoritic minerals from dust without 
loss of the carried xenon isotopes is too implausible; so this new idea was not sup- 
ported for isotopes of Xe. But the idea of such supernova dust had been firmly im- 
planted among the concepts of astronomy with radioactivity; moreover, the future 
may yet find evidence of fossil xenon. Several other extinct radioactive nuclei were 
predicted to exist within supernova dust (Clayton 1975a b), and many of these are 
now demonstrated to exist. 

Secondly, the 26 Al extinct radioactivity was discovered in meteorites. But its 
mean lifetime is too short (T26= 1.05xl0 6 yr) for any to survive the sequestering 
that xenology had suggested. This led to the second big change. The seeming im- 
possibility of maintaining the large measured quantity of interstellar 26 Al nuclei by 
nucleosynthesis prompted the idea that a supernova within the molecular cloud from 
which the sun was born had injected the solar 26 Al radioactive nuclei into the solar 
collapse just prior to the solar birth. Such a unique association would render the so- 
lar cl oud atypical of the ISM at large ( |Cameron and Truran 1977| Wasserburg et al 



1982). The substantial solar 26 Al abundance initially in meteorites became the stan 



dard bearer for a class of extinct radioactive nuclei that required special local pro- 
duction. Cameron and Truran ( 1977 1 suggested that far from being a coincidence, 
the supernova producing the solar 2b Al nuclei also triggered the collapse of the solar 
molecular cloud by the overpressure its shock wave brought to bear on that cloud. 
That model for live short-lived radioactivity in the solar cloud came to be referred 
to as the supernova trigger. But Cameron & Truran went too far. They supposed 
that the xenology -producing radioactivities ( 129 I and 244 Pu) were also injected by 
that supernova trigger. They ignored that the ISM already was expected to contain 
a hundred times too much 129 Xe for Reynolds' measurement, not too little, and ten 
times too much 244 Pu. So injecting those radioactivities seems to not be correct 
either. Nor does the idea of fossil 26 Mg nuclei from Al-rich interstellar dust rather 
than live 26 Al in meteorites seem correct. But the puzzles of extinct radioactivity had 
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come more sharply in focus, and these new ideas were instrumental for astronomy 
with radioactivity. 

Thirdly, additional extinct radioactive nuclei were discovered. So diverse were 
their lifetimes that the idea of exponential decay from a starting abundance became 
increasingly untenable. Then another new idea of great importance for astronomy 
with radioactivity appeared. Clayton ( 1983 ) pointed out that the concept of a homo- 
geneous ISM was at fault. He argued that supernova ejecta enter a hot phase of an 
ISM in which matter cycles on average among three phases, but that the sun was 
born from the cold molecular cloud phase. This distinction greatly modifies the ex- 
pectation of exponential decay in the ISM and of extinct radioactivities in molecular 
clouds. puss et al.| (|2009 ) provide a modern list of the many extinct radioactivity for 
which evidence exists. 

The revisiting of xenology and these three new developments left extinct radioac- 
tivity research changed conceptually. 



Neon-E 

A different noble gas, neon, provided an early but poorly understood evidence of 
Stardust in the meteorites. As with xenology, the 20 Ne/ 22 Ne isotopic ratio can be 
studied by incremental heating of a meteoritic sample in a mass spectrometer. David 
Black discovered that in the Ivuna meteorite the measured isotopic ratio 20 Ne/ 22 Ne 
decreased precipitously when the sample heating reached 1000°C, dropping to a ra- 
tio near three, much smaller than the common ratio near ten. Variation by a factor 
three in a trapped noble-gas isotopic component was unprecedented. This suggested 
that some unknown mine ral releases its neon gas near 1000°C, and that that mineral 



contains 22 Ne-rich neon. |Black et al.| dl972| suggested that this small 2U Ne/ Z7 Ne 



20\w22i> 



isotopic ratio was so bizarre that it can not be accounted for by conventional means 
within an initially homogeneous, gaseous solar system. Black suggested that inter- 
stellar grains that had formed somewhere where the 20 Ne/ 22 Ne isotopic ratio was 
smaller than three had survived the origin of the solar system and imprinted the 
1000°C temperature outgassing with its isotopic signature. He called this exotic 
neon component Ne-E. This was a radical far-reaching conclusion, perhaps the first 
of its kind based on sound analyzed data rather than on pure speculation. 

But where might this dust have formed? Arnould and Beelen] ( 1974} remarked 



that neon gas ejected from explosive He-burning shells of massive stars could re- 
semble Ne-E, but they offered no suggestion for forming carriers of Ne-E there or 
of otherwise getting such gas into the early solar system. Clayton] ( |1975b[ ) offered 
a more concrete proposal, arguing that supernova dust would condense before the 
ejecta could mix with circumstellar matter, roughly within the first year, and that 
22 Na (2.6 yr) produced by the explosion would condense as the element sodium and 
decay to daughter 22 Ne only after the grains had been grown. In that manner nearly 
isotopically pure 22 Ne could be carried within supernova dust into the forming solar 
system. This was an early prediction of the supernova-dust phenomenon in mete- 



orites. Clayton and Hoyle ( 1976 1 argued that nova dust provided another possible 
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source for Ne-E. Either radioactive venue would amount to a new technique for 
astronomy with radioactivity. 

The leading alternative to radiogenic 22 Ne is condensation of dust in a 22 Ne- 
rich red-giant stellar atmosphere. This can indeed occur, because during He burning 
the 14 N residue of the prior CN cycle can be converted to 22 Ne by two successive 
radiative alpha-particle captures ( |Arnould and Norgaard| [T978 ). The question arises 
whether this Ne is sufficiently 22 Ne-rich. And that required much more study of 
meteoritic specimens. But the possibility was experimentally strengthened by the 
discovery ( Srinivasan et al. 1978 ) that at least one Ne-E component was associated 
with i-process xenon from red-giant atmospheres. 

These were exciting developments in isotopic astronomy and for astronomy with 
radioactivity. This fresh new set of ideas intensified interrelationships among inter- 
stellar Stardust, the origin of the solar system, and meteorites. 



26 Al: Fossil or Injected Fuel 



During the mid 1970s it was discovered ( Gray et al. 1974| Lee et al. 1977[ ) that very 



refractory Al-rich minerals that apparently were among the first solids to form in the 
early solar system contained elevated isotopic ratios 26 Mg/ 27 Mg. Furthermore, the 
number of excess 26 Mg atoms within a mineral was shown to be proportional to 
the number of Al atoms. That correlation could be reproduced ( |Lee et aL~] |1977| > if 
the Mg had been initially isotopically normal and the initial Al contained a small 
component of radioactive 26 Al, whose later decay produced the excess 26 Mg atoms. 
The simplest way for this to have occurred is that the molecular cloud from which 
the sun had formed carried within it the large isotopic ratio 26 Al/ 27 Al=5x 10~ 5 at 
the time the planetary system began to form. The reason for calling that ratio large 
is that it far exceeded the ratio that had been anticipated to exist within the ISM. 
Had that ratio been applicable to the ISM as a whole, it would have corresponded 
to 29 M e of radioactive 26 Al nuclei - much more than was discovered in 1982 by 
Mahoney et al. by their detection of the 1.809 MeV gamma ray line that is emitted 
following its decay. 

The immediate issue was whether these radioactive 26 Al nuclei were actually 
alive in the early solar system or whether 26 Al was alive only when Al-rich dust 
formed in earlier galactic supernovae. If the latter option were true, the excess 26 Mg 
nuclei were fossils of radioactive 26 Al decay within interstellar Al-rich d ust. Exis 



tence of such 26 Mg-rich and Al-rich interstellar dust had been predicted ( Clayton 
|1975a|b[ ) in the first works to propose that dust containing evidence of extinct ra- 
dioactivity should condense within supernova ejecta. This was an exciting new idea 



in the astronomy of radioactivity. |Clay ton ( 1977a|b| l presented a solid-carrier model 



for excess 26 Mg nuclei in Al-rich minerals based on the fossil picture. In the first 
case, that of live 26 Al nuclei during the formation of the minerals, some explanation 
was required for why so much radioactive 26 Al should have been present in the so- 
lar matter when it could not have been a general property of the ISM. In the second 
case, fossil excess 26 Mg nuclei, why did interstellar Al dust containing abundant 
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radioactive 26 Al nuclei form within supernova ejecta and how did it participate in 
the chemical growth of the Al-rich minerals found in meteorites. In the first case, 
the heat of the radioactive decay of radioactive 26 Al nuclei within those planetesi- 
mals that had been the parent bodies for the meteorites would have been sufficient 
to melt large parent bodies, allowing them to differentiate chemically (as the earth 
has done). In that case 26 Al also played a role as fuel for melting of the meteorite 
parent bodies. In the second case, the excess 26 Mg nuclei represented a fossil of in- 
terstellar decay. Tension between these two possibilities was reflected in the phrase 
"Fossil or Fuel" within the title of the paper by Lee et al, who argued in favor of 
the live 26 Al nuclei in the solar gas. This debate raged for a few years but lost steam 
as chemical arguments for the growth of the mineral phases seemed increasingly 
likely to require live 26 Al nuclei. 

The struggle over the interpretation very short-lived extinct radioactivity was 
heightened later by the discovery of extinct 41 Ca in meteorites. Its abundan ce rel- 
ative to that of 41 K is only 41 Ca/ 41 K = 1.5xl0~ 8 ( jSrinivasan et al.[|l996| ); how- 
ever, even that small abundance looms large because its mean lifetime is but 741= 
0.144xl0 6 yr. Surely none can survive galactic nucleosynthesis to the time the solar 
system formed. So the supernova trigger injection would be needed to account for 
it as well as 26 Al. But fossil effects are also definitely possible. Huge 41 K excesses 
in Ca-rich supernova dust had been predicted two decades earlier (Clayton 1975b[ ). 

The history of documentation of extinct radioactive nuclei that had been alive in 
the early solar system was led for decades by G.J. Wasserburg and his laboratory 
at Caltech (jokingly called "the lunatic asylum" owing to their research on lunar 
samples). His laboratory was of such renown for its measurements using traditional 
mass spectrometry with extreme sensitivity that Wasserburg was chosen recipient 
of the 1986 Crafoord Prize for geosciences. Similar mass spectrometric research 
of great consequence for this problem was conducted by Giinter Lugmair in his 
laboratory at UCSD. 



Extinct Radioactivity and Mixing of Interstellar Phases 

Xenology soon gave way to a host of newly discovered extinct radioactive nuclei. 
129 I and 244 Pu could seen to be members of a larger class of now extinct parents. 
A very special class contains those nuclei having mean lifetimes long enough that 
some of that radioactive abundance might survive from galactic nucleosynthesis 
but short enough that their survival live to the time of the solar system depends on 
the timing of interstellar mixing. That class of nuclei also includes 129 I and 244 Pu. 
Their mean lifetimes are, in ascending order: f( 60 Fe)= 2.1(Qx 10 6 yr; T( 53 Mn)= 
5.34xl0 6 yr; r( 107 Pd)= 9.38xl0 6 yr; T( 182 Hf)= 13.8xl0 6 yr; t( 129 I)= 23.5x10 6 
yr; t( 244 Pu)= 115x10 6 yr; and r( 146 Sm)= 149xl0 6 yr. The relative abundances 
that they possess and the three-isotope correlation plots by which they are measured 
are discussed in Chapter 6. What became abundantly clear is that simple exponential 



The mean lifeof 60 Fe is now known as 3.8 My 
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decay from a set of starting abundances is not consistent with the data. Something 
more complex is at work. 

Four layers of interpretation The new data and new creative ideas led to a 
picture of cosmic radioactivity that is interpreted on four layers of complexity. Like 
the proverbial Russian doll, one opens each layer only to reveal a more complex 
version inside. What are these four layers? 



1. Mean Galactic Nucleosynthesis and the well-mixed ISM: Using Eq. |2.21| for the 
mean interstellar ratio of a radioactive abundance to that of a reference isotope 
gives that mean abundance ratio to be 

N( A Z) /N( A *Z) = P{A Z ) /P{ A *Z) (k+l) T Z A© (2.22) 

If the halflife of species A Z is long enough that the ISM can be considered to 



be well-mixed, it may be interpretable by the mean expectation of Eq. 2.22 An 
example might be the uranium isotope ratio, since both have halflives of 10 9 yr 
or greater. The ISM probably mixes well during 10 9 yr. 
2. Mean Galactic Nucleosynthesis and the ISM Phases: If the lifetime Tz is of or- 
der 10 8 yr or less, the abundances of that radioactive nucleus may differ in the 



different phases of the ISM. The three-phase mixing model (Clayton 1983 i that 
describes this in the mean is described below. The idea is that supernova ejecta 
appear initially in the hot ISM, and considerable mean time is required for it 
to work its way via mass exchange between phases into the molecular clouds 
where the sun might have formed. That delay reduces radioactive abundances, 
depending upon their lifetimes. This idea regards each phase of the ISM as hav- 
ing the mean concentration of radioactive nuclei everywhere in that phase; but 
the three phases have differing mean concentrations of the radioactivity, with the 
molecular-cloud abundance being the least. This picture divides the mean ratio 



Eq. 2.21 into differing steady-state ratios in each of the three phases. An example 
might be the iodine isotope ratio, since its halflife of 17 My is shorter than the 
expected phase change times (perhaps 50 Ma) for the ISM phases. Thus each 
phase will in the mean have different 129 l/ 127 l ratios. 

Mean Galactic Nucleosynthesis and Deviations from mean ISM mixing times: 
All galactic samples of a given ISM phase will surely not be identical. Differ- 
ences occur because nucleosynthesis does not follow the galactic mean rate at 
all locations, or because the interphase mixing times may vary from one place 
to another. Separate portions of a phase will scatter about that mean for that 
phase, perhaps by sizeable amounts. Astronomical observations of very old, very 
metal-poor stars reveal this phenomenon clearly (Burris et al. 2000). The iodine 



isotopes can illustrate how this works. Consider the molecular cloud from which 
the sun formed. It may have existed in a portion of the galaxy whose hot phase 
received less than the average amount of new nucleosynthesis during the past 
24 My (the 129 I mean lifetime). Thus the 129 I concentration locally would be less 
than its galactic mean, whereas stable 127 I would have the same concentration 
everywhere. Furthermore, the time delay for the hot phase to be admixed into 
the solar molecular cloud may locally have exceeded the galactic mean delay, in 
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which case the 129 I concentration in the solar molecular cloud would have de- 
cayed to a smaller value than typical of that phase. This layer of interpretation 
deals with specific local galactic workings that differ from the galactic mean. 
4. Nucleosynthesis from a single neighboring supernova: The solar concentration 
of a short-lived radioactive nucleus may have been the result of a single galactic 
supernova. The shorter the halfiife, the more likely it becomes that a neighbor- 
ing dying star is needed to create it. In such a case the concept of mean galactic 
nucleosynthesis rate has no relevance. Short-lived 26 Al is historically the prime 
example of this situation. Variations on the supernova injection model then ap- 
pear: When did that supernova occur? With what efficiency did its ejecta admix 
into the solar molecular cloud? How long was required for that mixing to occur? 
Such questions are so specific that they are frustrating. Almost any answer seems 
possible, depending upon hydrodynamic modeling for its credibility. 

Finally, some radioactive nuclei may have multiple causes for their solar abun- 
dances. Each layer of interpretation may contribute. Take 244 Pu as a likely example 
(t( 244 Pu) =115 Ma). A portion of the solar 244 Pu abundance may survive at the level 
expected from mean galactic nucleosynthesis (layer 1). Layer 2 may have reduced 
that somewhat, depending on the mean interphase mixing times. Layer 3 describes 
a non-average 244 Pu solar concentration owing to atypical (non-mean) local nucle- 
osynthesis rate and mixing parameters. Finally, a portion of solar 244 Pu abundance 
may have resulted from a single nearby supernova at the time of solar birth. It will 
be clear that scientific understanding must seek the best picture for simultaneously 
fitting all of the radioactivities. Each individual abundance will be inadequate by 
itself to determine the physics of its solar presence. But a model that fits them all is 
a serious contender for the truth. To advance more into this topic, consider the mean 
interphase mixing model. 

Mean mixing model for ISM radioactivity: Astrophysics provides this novel 
example of the layer-2 mathematics of radioactive decay. The distinct ISM phases 
have vastly differing temperatures and densities. These phases mix with one an- 
other on timescales that are longer than the shortest of the extinct radioactivities 
but shorter than the longest. As a consequence, the concentration X of a radioac- 
tive nucleus will differ from phase to phase, with the differences greater for shorter 
halfiife. Solar-system samples revealing the presence of once-live radioactive nu- 
clei are obtained from meteorites and other planetary objects formed in the solar 
accretion disk. Measurements of solar samples therefore are of the concentrations 
in the cold molecular cloud (one of the phases) from which the sun was born. The 
key point is that radioactive concentrations in that molecular-cloud phase should be 
smaller than in the mean ISM, which includes all of its phases. This consequence of 
the mixing times among the phases, taking into account their levels of radioactivity, 
was devised by Clayton (1983). 



A simplified two-phase model To focus on mathematical aspects of radioac- 
tivity rather than on aspects of astrophysics, consider first a simplified model that 
illustrates the essence of the problem. Suppose that the ISM consists of but two 
phases. Stars are born in phase 1, but the freshly synthesized radioactivity is ejected 
from supernovae into phase 2. Let those two phases have equal mass M and ex- 
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change matter with the other phase on a timescale T. The mass exchange rate is 
then (dM/dt) ac f, = M/T. Matter from phase 1 joins phase 2 at that rate, and con- 
versely matter from phase 2 joins phase 1 at that same rate. The masses of each 
phase then remain constant and equal. 

In reality, spatial inhomogeneity also exists, depending on where in phase 2 
the fresh radioactivity is created and how long is required for its homogenization. 
Roughly 30,000 supernovae occur per My in the Milky Way, and they have simi- 
lar nucleosynthesis yields, so an almost homogeneous injection phase occurs faster 
than the output of a single supernova can be spread uniformly. Even so, concentra- 
tion differences and isotopic differences will exist spatially in that injection phase. 
The best that can be done without a specific calculation defined for a specific config- 
uration is to calculate the average concentration in the injection phase. To this end 
consider the injection phase to be uniformly mixed at all times (Layer 2). With that 
assumption the average difference in radioactive concentration between the phases 
can be calculated with a simple mathematical model. The answer can not be as- 
sumed to have applied exactly to the unique case of solar formation. Nonetheless, 
such a calculation reveals what is anticipated on average, without inhomogeneous 
effects. 

Returning to the simplified two-phase model for those average concentrations, 
the mixing rate between them is (dM /dt) exc i, = M/T. Then with X\ being the con- 
centration (grams of X per gram of ISM) of a radioactive nucleus in phase 1, the 
star-forming phase, its total radioactive mass MX\ changes as follows: d/dt(MX\) = 
MdX\/dt +X\dM jdt = MdX\jdt because the second equality is a consequence of 
the unchanging mass M of each phase. How shall one evaluate dX\/dtl It changes 
owing to two effects: first, the radioactive decay of X\\ second, the mass exchange 
with phase 2 having concentration X^. One writes for those two terms MdX\/dt = 
—MXi/t+ (Xz —X\)M/T where the first term is the radioactive decay rate with 
mean lifetime T. But the left-hand side must also vanish because dX\ jdt = in the 
steady state. This reasoning yields at once X\ = t/(T + t)X2 This factor expressing 
the ratio of the concentrations can differ significantly from unity. Suppose the mix- 
ing time T=100 Ma. Then for 129 L X\ = (23.5Ma/123.5Ma)X 2 = 0.19X 2 . The 129 I 
concentration in star-forming clouds is then fivefold smaller than in the injection 
phase. This reduction is even greater for a nucleus having smaller mean lifetime. 
But for very long lifetimes, l/(T + t) approaches unity, both phases having the 
same concentration. 

Both X2 and X\ can be individually evaluated from the requirement that their av- 
erage (since the two phases have the same mass) must equal the mean expected ISM 
concentration for that nucleus (Eqs. |2.20| and |2.2l) . Chapter 6 includes a table of the 
known extinct radioactive nuclei and their abundances in the early solar system. 

The significance for differing ratios of extinct nuclei between the two phases is 
that those differences do not reflect the simple expectation of exponential decay. The 



ratio t/(T + t) is very far from exponential in dependence on lifetime T. Clayton 



( 1983 ) introduced that effect when comparing the solar abundance levels of differing 
radioactive nuclei. This simplified two-phase mean model illustrates that idea. 
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Mean Life (Years) 



Fig. 2.2 Calculated surviving levels of the extinct radioactive nuclei in the solar cloud at the time 
of solar birth are expressed as a ratio to a neighboring stable isotope of the same element. Galactic 
steady state modefied by the continuous interphase mixing is shown for several different mixing 
times, (from Huss et al. 2009 1 . 



Clayton three-phase model The initial (Clayton 1983) introduction of this 



physical idea for interpreting relative abundances of extinct radioactive nuclei in this 
manner actually suggested a three-phase ISM rather than only two. Those masses 
and phases were: 

(1) Mi is the molecular-cloud mass in which stars form. It exchanges mass with 
M 2 . 

(2) M2 is the mass of large HI complexes that can not be disrupted by the super- 
nova shock waves that frequently traverse the ISM. Phase 2 exchanges mass with 
M\ with lifetime Ti and also with M3 with lifetime Ti- 

(3) M3 is the mass average HI clouds that are sufficiently small to be disrupted 
by passing supernova shock waves and which therefore are part of this warm neutral 
phase. Phase 3 exchanges mass with phase 2 with lifetime Ti- 
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It is not clear that this three-phase model is more realistic than the simpler 
two-phase model; but the need for fresh radioactivity to mix from its injection 
phase through a second phase before admixing with molecular clouds does amplify 
the distinction between radioactive lifetimes. Future progress with this multiphase 
idea, which definitely is significant for the extinct radioactivity problem, will re- 
quire more modern astrophysical work on the differing phases, their masses, and 
how effectively they exchange mass with one another. Finally, it may be that time- 
dependence must be introduced. The growth of each phase, rather than static mean 
masses for each, needs evaluation. Over long times the global masses of each phase 
will be static; but their cyclic growth and disruption may be an essential part of the 
problem. 

An insightful paper in regard to these issues has recently appeared ( |Huss et aL| 
2009 ). Those authors present current data for the observed abundances of each ex- 
tinct radioactivity. They also review the sites of nucleosynthesis of each nucleus. 
Then they evaluate the mean abundance expected within the standard model of 
galactic chemical evolution moderated by the interphase mixing described above. 
Fig. 2.2 (from [Huss et al. , 2009 ) displays the expectation as a function of the mean 
decay lifetime for a wide range of interphase-mixing timescales and compares those 
results to those for two differing free-decay intervals (the solid curves). Note that 
interphase mixing fits more extinct radioactivities than does either free-decay inter- 
val;but the four short-lived isotopes are too abundant for any galactic scheme, ap- 
parently demanding injection from an adjacent supernova. Keep in mind that what 
is being examined by these mathematical procedures is a fundamental property of 
interstellar radioactivity, couched for convenience and transparency in terms of re- 
alistic analytic representations of the true galactic mixing processes. 



2.2.5 Stardust: Radioactivity in Solid Samples ofPresolar Stars 

A presolar Stardust grain is composed only of atoms from a single donor star. The 
word Stardust applies to only those grains. It is used as a scientific word rather than 
a poetic word. It may be hard to accept that the isotopic composition of elements in 
an interstellar grain attests to it being a solid chunk of a single star. An average ISM 
grain has suffered a lengthy residence time within the interstellar medium before 
being incorporated into a new stellar accretion disk. Do not high energy interstellar 
ions implant within the grains as well as sputter them? Can interstellar chemistry 
add more atoms to mineralized grains? Would not grain ejection from a protostellar 
disk after chemically adding atoms to it within the hot dense disk structure also 
produce mineralized grains? Many questions flood the skeptical mind. Such doubts 
are reasonable considering the incomplete knowledge of these and other processes. 

Stardust grains are not average interstellar grains. Those grains themselves pro- 
vide the best answers to the doubts expressed above. Their well ordered crystal 
structures suggest thermal condensation at high temperatures, as would be expected 
within hot but slowly cooling gases leaving a star. Annealed crystals would not, 
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in contrast, be expected from accreted interstellar atoms. The dramatic isotopic ra- 
tios measured within individual Stardust grains, obtained from nearly homogeneous 
grains having large numbers of atoms (10 6 to 10 12 , say), strongly suggests that the 
condensation was from a gas having the extremely non-solar isotopic composition of 
the Stardust grain. The SiC grains, one of the most abundant of all presolar Stardust 
grains, were assembled from C atoms having observed isotopic ratios in the grains 
between 12 C/ 13 C=3 to 5000, whereas hot chemistry in an accretion disk would shift 
ratios toward interstellar norms, 12 C/ 13 C=89 in the case of the solar accretion disk. 
Similar evidence appears in isotopic compositions of N and Si, to name only the 
most studied of the elements. Those ratios attest to an extreme condensation envi- 
ronment wherein isotopic composition is not that found in the ISM. These isotopic 
studies of individual Stardust grains have been made using SIMS (secondary-ion 
mass spectrometry) in terrestrial laboratories (Bernatowicz et al. 1997 1. 

What followed from the discovery of Stardust was nothing less than a revolution 
in astronomy. Solid pieces of stars are now characterized routinely in terrestrial labo- 
ratories, including isotopic analysis of many elements with a precision not attainable 
at telescopes. Distinguishing isotopes in stellar spectra is very difficult. Three scien- 
tists dominated the early development of this field experimentally. Robert M. Walker 
committed his laboratory at Washington University to development of new technol- 
ogy capable of better laboratory study of primitive solar system samples. He fore- 
saw the capability of the sputtering ion probe, and recruited Ernst Zinner to devote 
a decade to improving and studying that technique. Their laboratory at Washington 
University (St. Louis) was then ready to study the newly isolated meteoritic grains, 
which they documented to be Stardust. Zinner and his many students led and domi- 
nated the new astronomy that blossomed from their isotopic analysis. An engaging 



history of this has been written in tribute to Walker and Zinner by McKeegan et al. 



(2007 ). Thirdly, Edward Anders inspired his associates at University of Chicago to 
the detective-story pursuit that isolated Stardust from the bulk meteorite rock. They 
termed their approach burning down the haystack to find the needle. It tracked the 
chemical carriers within the meteorites of isotopically anomalous neon and xenon 
through a sequential dissolving the bulk of the rock in strong acids ( |Amari et al.| 
1994} and finding that the undissolved residue became increasingly anomalous. The 



carriers of anomalous noble gas were, fortunately, carbonaceous, and therefore did 
not dissolve in acid. From these residues they were able to extract the individual 
Stardust grains. This too was a profound adventure in scientific exploration. Finally 
it may be noted that predictions of the existence of isotopically anomalous Stardust 
had existed for a decade prior to their experimental discovery. 

Within oxide grains ( |Nittler et al.| |1997| l stunningly abnormal families of cor- 
related O isotopes rule out any growth within mixed ISM. Equally stunning was 
the almost pure isolated s-process isotopic compositions that were predicted for 
noble xenon in red giant Stardust (Clayton and Ward |1978| l and discovered (Srini- 



|vasan et"aL| |1978| l in bulk carbon-rich residues of acid dissolution of meteorites. 
These Xe s -rich residues are dominated by collections of SiC grains. By contrast, in- 
terstellar bombardment would implant normal Xe isotopic composition. Neon iso- 
topes also revealed an exciting story. Individual mainstream SiC grains have been 
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shown to carry almost pure ^-process Mo, a quite reactive trace element rather than 
a noble gas, showing that the Mo-containing condensing gases had not mixed with 
the ISM prior to condensation. Isotopic evidence of Stardust abounds. By careful 
consideration of the entirety of the known properties of Stardust grains, researchers 
become confident that they are indeed solid samples of stars that died before the sun 
was born. 



Stardust predicted 

With the realization that computer models of massive-star supernovae generate 
shells of abundant intermediate-mass elements, Hoy le et al.| ( |1970| l argued that the 
adiabatic expansion and cooling of these newly synthesized nuclei within the super- 
nova interior should be accompanied by condensation of a potentially large amount 
of silicate dust. The ejected gas cools to below 2000K after only several months 
when the density is still quite high. Their goal had been to account for the interstel- 
lar dust mass of silicates and of graphite. It did not occur to them that such grains 
could be identified by their unusual isotopic compositions. Nor did they suggest that 
they might be found within meteorites. They simply tried to explain the observed 
existence of dust in astronomy. 

Clayton (1975a]b 1978] ) advanced exciting observational considerations on the 



grains' condensed radioactive nuclei. It was suggested that grains from the supernova- 
condensation process would be identifiable by their extreme isotopic signatures, in- 
herited from the isotopic compositions of those supernova shells. Radioactivity is 
prominent within supernova shells. Short-lived radioactivity was proposed to de- 
cay within each grain during its interstellar residence. So the ISM should contain 
everywhere numerous interstellar grains from the very large number of presolar su- 
pernovae. Also proposed was a chemical memory model (Cla yton| |1977a|b( [1978 , 
1982) > in which supernova grains were incorporated during later growth in the solar 



accretion disk of the larger solar-system solids. Many of those solids were incor- 
porated into meteorites. That picture suggested potential explanations for several 
isotopic anomalies that had been discovered in meteorite solids; but it was rejected 
initially by meteoriticists, who focussed their doubts on the potential explanations 
of xenon anomalies from extinct radioactivity. Clayton and Wickramasinghe ( 1976} 



and Clayton| ( 1979| l; Kozasa et al. ( 1989 1 demonstrated that a sufficient number of 



collisions of supernova-interior atoms will occur with any grain nucleation during 
expansion of the supernova interiors to grow grains greater than one micrometer in 
radius. None of these works was so optimistic as to predict that individual supernova 
grains would be found intact within meteorites; but they did establish the ideas by 
which those discoveries could be recognized. Radioactivity lay at the heart of these 
predictions. 

In a paper outlining a system for the different types of precondensed matter in 
the early solar system, the supernova condensates were named SUNOCONs and 
the term STARDUST was restricted to hot thermal condensation during mass loss 
from other stars (Clayton 1978| l, primarily red giants. NEBCON was suggested for 
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grain mantles grown by nebular sticking of atoms and molecules to preexisting dust. 
These names have not found the favor of usage, despite describing a theoretical pic- 
ture of what to expect in the early solar system. Only the term Stardust is commonly 
used, and applied to all types of high-temperature thermal condensates from stars. 
The title Precondensed matter: key to the early solar system of Clayton ( 1978 ) ex- 
plicitly contradicted that epoch's favored picture in which the initial solar system 
was initially hot and totally gaseous and condensed its solids from that solar gas. 

The predictions of isotopically extreme Stardust were, after some initial contro- 
versy, wonderfully confirmed by the discovery of supernova Stardust and the later 
documentation of chemical memory within large solids, earning broad acceptance. 



Stardust discovered: Fossil Extinct Radioactivity 



Secondary-ion mass spectrometry (SIMS) of ions sputtered from isolated single me- 
teoritic grains by a focused ion beam identified Stardust in the late 1980s. A history 
may be read in Anders a nd Zinner| ( | 1 99 3 | l . First to be discovered ( |Bernatowicz et al.] 
1987 IZinneret al.||1989) was silicon carbide (SiC), with isotopic patterns that iden- 



tified it as thermal condensation in matter flowing away from the photospheres of 
asymptotic-giant-branch red giants. This sensational identification of Stardust is now 
beyond doubt, buttressed by a huge number of experimental and theoretical studies. 
These carbide grains can condense only after the composition changes caused by the 
third dredgeup of carbon, which eventually cause carbon to be more abundant than 
oxygen in the envelope. As long as O is more abundant than C, the CO molecule 
oxidizes (combusts) all carbon. Those SiC grains condensed after the carbon-star 
transition contained 26 A1/ 27 A1 near 10~ 3 based on the excess 26 Mg atoms within 
their magnesium-bearing phases. When supernova SiC X grains were later identi- 
fied as a subfamily of SiC Stardust, many revealed 26 A1/ 27 A1 near 0.1 when they 
condensed within the expanding supernova interior. 

Oxide St ardust too was iden tified and classified into isotopic families a few 
years later (Nittler et al. 1997 1. Its most common pattern is n O-richness owing 
to secondary production in the parent star of 17 from initial ie O during hydrogen 
burning . At about the same time supernova Stardust was identified unambiguously 
( |Hoppe et aL] |1996[ |Nittler et al] |1996[ ) using the prediction ( |Clayton| |1975b| l that 



excess^Ca would exist in the Ti-bearing phases within grains owing to the conden- 
sation of radioactive 44 Ti within a year after the explosion. Because the elemental 
ratio Ti/Ca is for crystal-lattice reasons larger than solar in SiC grains, the supernova 
SiC grains possess very large isotopic excess at 44 Ca. These were found in the fam- 
ily of SiC grains called X grains, which were already suspected of being supernova 
grains on the basis of deficiencies in the heavy iso topes of Si and C and on excesses 



of 15 N (Alexander et al 



1990 



Amari et al. 



1992) 



The extremely large isotopic excesses of 26 Mg and 44 Ca in these SiC X grains 
were the most dramatic discoveries of fossil extinct radioactivity in Stardust. They 
demonstrated that radioactive 26 Al and 44 Ti had been quite abundant when the grains 
condensed during the supernova expansion. Predictions had also suggested that fos- 
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sil abundances in interstellar grains might contribute to excesses of the daughter iso- 



tope in solids grown later from interstellar grains (chemical memory). But Lee et al. 
( |1977| | had emphasized grounds for believing that the excess Mg within solids 
made in the solar system reflected instead live 26 Al at the time the solids were made. 
Both fossil and live 26 Al now appear to have existed. These discoveries dramatically 
spotlighted rich new roles for the astronomy of radioactivity in early solar system 
chemistry. 

Many other fossil extinct radioactive nuclei have been added to the observed list. 
In supernova Stardust these include excess 22 Ne owing to 22 Na decay, excess 41 K 
owing to 41 Ca decay, and excess 49 Ti owing to 49 V decay within the grains (see 
|Clayton and Ni tiler 2004). All had been predicted to exist in supernova dust, but 
their discovery surprised and delighted isotopic chemists and astronomers alike. 

Following the exciting discovery of presolar Stardust, its study has emerged as a 
new area of astronomy. In just 23 years its existence has evolved from a bewildering 
new discovery into several new techniques for measuring isotopic abundance ratios 
with high precision in stars. This has been especially important for astronomy with 
radioactivity, because in Stardust the level of extinct radioactivity reflects its abun- 
dance during nucleosynthesis in stars, whereas in solar system samples it reflects 
the level of its survival. Precisely measured isotopic ratios for four to eight chem- 
ical elements endows each gemlike refractory mineral with significance for stellar 
structure and evolution, and for the chemical evolution of the Milky Way during the 
epoch 7-5 Gyr ago, and for new insight into nucleosynthesis ( |Clayton and Nittler] 



2004). This rich harvest is compromised only by the fact that the donor stars can not 



be seen, because they died before our solar system began. Their Stardust bears no 
label, save that of their measurable properties. The nature and evolutionary state of 
the donor stars must be ascertained from the detailed properties, primarily isotopic, 
of each grain and of the way each grain fits into the spectrum of the thousands of 
Stardust grains that have been analyzed to date. The rapidly growing numbers of 
analyzed grains allows evolutionary trends within their distinct families to be iden- 
tified. 



Radioactivity and Chemistry of Condensation in Supernovae 

An unanticipated aspect of astronomy with radioactivity lay in the role of radioac- 
tivity in the chemistry of the condensation process. The supernova interior offers 
a unique laboratory for condensation physics. It guarantees that chemistry begins 
with gaseous atoms, with no trace of previous molecules or grains. Gamma rays 
and their Compton-scattered electrons bathe the supernova core. By disrupting the 
CO molecule, they cause its abundance to be very much smaller than expected from 
states of chemical thermal equilibrium at the ambient temperature. The small CO 



4 Readers wishing familiarity with these topics can best consult later chapters of this book, or 
|Clayton and Nittler (2 004), a nd Astrophysical Implications of the Laboratory Study of Presolar 
Materials |Bernatowicz et al. 1997); the Handbook of Isotopes in the Cosmos by |Clayton| l |2003) 
also presents many astrophysical consequences. 
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molecular abundance enables other non-equilibrium paths to the condensation of 
carbon-bearing solids. This disequilibrium can be considered to be another aspect 
of radiogenic luminosity causing a higher degree of excitation in the gas than would 
be expected in the absence of radioactivity. 

Because the supernova core is hydrogen-free, chemical condensation routes uti- 
lizing H are not relevant. Polycyclic aromatic hydrocarbons, for example, do not 
come into play save in the envelope; but 90-95% of the mass of ejected Mg, Al and 
Si emerges in the H-free core. And yet abundant dust condensation is observed to 
have occurred in SN 1987A (|Wooden et al) [1993] |Colgan et aL ] [T994j), in Cas A 



(|Arendt et al. 1999 Dunne et al. 2003| l and in the Kepler remnant ( Dunne et al 



|2003| l. Apparently a few solar masses condensed in Cas A and about one solar mass 
in Kepler, requiring high condensation efficiency for Mg, Al and Si and even for 
carbon and thereby establishing supernovae as a major contributor to the budget of 
thermally condensed interstellar grain cores. 

The traditional guideline to condensation of solids had been to assume that chem- 
ical equilibrium applies during the expansion and temperature decline of the super- 
nova interior. Formulae yield the equilibrium condensed masses of differing min- 



erals (Lattimer et al. 1978 Ebel et al. 2000 1. This approach can not yield grain 



size but only the total condensed masses. The assumption of thermal equilibrium in 
the solid phase exaggerates the ability of chemical reactions to maintain chemical 
equilibrium within the solids as the gas cools. Because of the rapid fall of density 
and temperature during the expansion of the supernova interior, thermal condensa- 
tion must be accomplished within about two years, too fast to maintain equilibrium. 
A nonequilibrium theory of condensation based on nucleation theory followed by 



subsequent growth has been developed by Kozasa and collaborators ( Kozasa et al. 



[1987] [1989] [1991] |Todini et al.] [2001) . Their method identifies a key molecular 
species whose abundance controls the condensation. However, their questionable 
assumptions concerning gaseous mixing at the molecular level and on an outdated 
role for the CO molecule render their results of questionable value. 

The nucleosynthesis problems posed by isotopic ratios within individual super- 
nova grains can not be decoupled from physical questions about the chemistry of 
their condensation. Condensation chemistry is an essential aspect of astronomy with 
radioactivity. Because no single supernova mass zone is able to satisfy the isotopic 
ratios found in supernova grains, it has long been clear ( |Amari etal] \991) that some 
type of mixing before condensation is needed to produce their chemical and iso- 
topic compositions; but it is not clear whether that mixing represents atom-scale 
gaseous mixing in the very young remnant or transport of a growing grain from one 
composition zone into another. But the biggest question is whether the requirement 
that the C abundance be greater than O abundance in order to condense SiC and 
graphite within supernovae is a valid requirement. For their discussion of supernova 
Stardust, for example, Travaglio et al. (1999) took the view that mixing occurs at 
the atomic level, prior to condensation, and that only C>0 portions of those mix- 
tures support condensation of carbon. They arbitrarily and instantaneously mixed 
regions having C<0 with the He shell having C>0 to find mixing fractions that 
retained C>0 overall. Kozasa et al. ( 1989) performed similar thought experiments. 
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Such treatment is deeply flawed, however, because it seems certain that gases can 



not mix at the molecular level within a few years time (see Deneault et al. 2006 



Fryxell et al. 1991 1. The so-called mixing calculated in hydrodynamic studies of su- 
pernovae, on the other hand, represents fluids of one composition exchanging places 
in the homologous expansion with overlying fluids of a different composition. This 
is a different use of the word mixing, and has introduced considerable confusion 
into published works. Mixing at the molecular level requires very much more time 
than the year that is available before condensation must occur. Later the expanding 
remnant is too dilute, its density too small, for condensation. 



To circumvent the equilibrium CO trap for carbon, Clayton et al. ( 1999} , Clay- 



|ton et al.] ( |2001| l, and |Deneault et al7| ([2006 ) called upon its dissociation by the fast 
compton electrons energized by supernova gamma rays. They advanced a kinetic 
theory of graphite growth and calculated its consequences in detail after advocat- 
ing a specific nucleation model. The key idea is that small graphite test particles, if 
placed in a hot gas of C and O atoms, will associate with free C atoms faster than 
they can be oxidized by free O atoms. This is a property of reaction cross sections 
rather than of chemical equilibrium. Even though oxidation of carbon would surely 
be the ultimate end given adequate time, the expansion will terminate the chemistry 
after about two years with large graphite grains remaining. The graphite essentially 
is a metastable state of carbon. This theory is supported by observing that super- 
nova 1987A ejected only 10~ 3 M & of CO molecules instead of the 0.1 M Q of CO 
molecules that is first formed by association reactions in the hot gas before radioac- 



tive disruption reduces its abundance (Dalgarno et al.| 1994 1995 Gearhart et al. 
[1999). 

Similar issues probably surround the condensation of supernova SiC. It seems 
plausible that radioactive liberation of free C atoms from CO molecules will also 
facilitate the condensation of SiC in O-rich gas; but this is hotly debated. Although 
a kinetic route to SiC condensation has not been layed out, Denea ult et al.| (|2"003b ) 
have formulated a physical description of the ejecta enabling them to make sev- 
eral relevant conclusions derived from assuming that the radioactive CO-disruption 
mechanism is the correct key to SiC growth. They present animations of a hydrody- 
namic calculation showing that a reverse shock wave launched toward inner zones 
by the radially increasing value of the product pr 3 in the H envelope compresses 
a dense shell near radial mass coordinate m = 3M , where Si and O are the most 
abundant elements; but some C abundance remains there for possible carbon chem- 
istry. They propose that SiC condenses in that zone, and they detailed several other 
physical processes that may produce the observed grain compositions. Mixing of 
a new type during condensation also occurs if the reverse shock generated by the 
radial expansion colliding with the presupernova wind arrives at the condensation 
zone at m = 3M Q between six months to a year after the explosion, because that 
reverse shock slows the gas and forces the partially condensed SiC grains to propel 
forward through the decelerating gas into more 29 :30 Si-rich regions, giving perhaps 
a new interpretation for mixing during condensation. 

Despite many uncertainties, it now appears certain that supernova grains studied 
by isotopic analysis will provide, through details of condensation chemistry, a new 
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sampling spectrum of young supernova interiors, just as have gamma-ray lines and 
hard X rays. If so, a very rich but complex discipline of astronomy with radioactivity 
will follow. 



Astronomy with Radioactivity Today 

This introductory chapter has focussed on the history and the key physical ideas 
of astronomy with radioactivity. Effective research depends on a clear grasp of an 
interdisciplinary set of its central ideas. History itself often provides the best exam- 
ple by which the physical idea can be first grasped. But the remainder of this book 
looks forward, not backward. Every aspect of astronomy with radioactivity today 
involves grappling with a host of technical details. This book attempts to bring the 
reader to that point. Each aspect of astronomy with radioactivity today also involves 
grappling with the entire world of astronomy. Astronomy itself is many disciplines, 
and not even astronomers are always comfortable outside their own astronomical 
technique. How much harder it is to place the complex manifestations of radioactiv- 
ity into a continuously changing astronomical fabric. But this is the direction of all 
scientific progress. 
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